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ABSTRACT 


Di;  ital  Simulation  of  Ship  Propulsion  Trains  Utilizing 
Gas  Turbine  and  Diesel  Prime  Movers 

Ernest  R.  Freeman 

A  computer  simulation  of  some  power  plants,  transmission 
systems,  and  propellers  has  been  developed. 

Theoretical  models  have  been  established  mainly  from 
published  models-  for  the  following  components. 

1.  prime  movers 

a.  aircraft  derivative  gas  turbines 

b.  low  speed  diesel  engines 

c.  medium  speed  diesel  engines 

2.  transmissions 

a.  direct  drive  (low-speed  diesel) 

b.  reduction  gear  with  and  without  reversing 
clutch  ar r an ge men t 

3.  propulsors 

a.  conventional  fixed  pitch  propellers 

b.  controllable  reversible  pitch  propellers 
The  computer  simulation  utilized  the  Dynamic  System 

Simulation  Program  (DYSYS)  on  an  Ihh  1130  computer.   DY3YS 

was  developed  by  the  Department  of  Mechanical  Engineering 

at  M.I.T.   Representative  accelerations  and  crashback 

maneuvers  were  computed  ising  the  various  propulsion  systems'. 

Several  critical  situations  were  established  involving 

the  sequencing  of  power  plant,  transmission,  and  propulsion 

transients. 

Thesis  Supervisor:  A.  Douglas  Can-ichael 
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INTRODUCTION 

The  rapid  advances  made  in  the  field  of  auto?atic  controls  since 
World  War  II  have  found  more  and  more  applications  in  the  area  of  ship 
propulsion.  One  of  the  principal  driving  forces  behind  this  is  the 
desire  to  operate  plants  safely  and  at  or  near  their  "optimal"  operating 
point*  The  safe  coordinated  operation  of  some  ship  power  plants  con- 
sisting of  high  performance  prime  movers,  sophisticated  couplings,  and 
propellers  is  no  trivial  matter*  This  is  especially  true  during  periods 
when  fast  ship  response  is  desired,  e.g.,  emergency  stopping  or  acceler- 
ating « 

To  be  able  to  control  a  propulsion  plant,  on*  must  be  able  to 
model  it  and  kn<;v  how  it  behaves  under  various  disturbances  or  changes 
in  inputs.  The  alternative  is  a  costly  experimental  program  for  each 
system.  The  problem  for  the  control  engineer  now  becomes  one  of  developing 
or  choosing  a  model  that  suitably  describes  the  characteristics  of  the 
particular  power  train  cociponents.  These  .component  models  may  be  either 
physical  or  functional.  A  physical  model  may  accurately  describe  a 
component's  characteristics  using,  for  example,  Newton's  Laws  of  Motion v 
thermodynamic  relationships,  conservation  of  mass  relations,  etc.  These 
provide  a  valuable  insight  of  the  physical  processes  occurring  in  the 
operation  of  a  component,  and  enable  one  to  perceive  what  simplifications 
can  be  made  to  a  model.  Some  components,  on  the  other  hand,  may  prove 
to  be  very  non-linear  and  to  describe  their  operational  characteristics 
using  physical  laws  may  result  in  mathematical  expressions  which  are  not 
suitable  for  use  in  a  simulation  study.  In  this  situation  a  simpler 
mathematical  expression  is  developed  which  closely  approximates  the 
component's  characteristics  over  a  reasonable  range  of  inputs. 


IMTR01XTCTJ0M  ( Cont  'd ) 

The  philosophy  behind  the  study  described  in  this  thesis  has  been 

to  keep  the  models  simple,  but  yet  reasonable*  To  quote  p.  J.  Kiviat, 

who  prepared  a  RAND  memorandum  on  digital  simulation  modeling  concepts, 

"the  model  should  only  be  as  detailed  as  is  necessary  to 
answer  the  question  at  hand;  it  should  be  so  designed, 
however,  so  it  can  be  expanded  to  include  more  detail 
without  inordinate  cost  in  those  model  areas  which  have 
a  high  probability  of  becoming  subsequent  subjects  of 
concern. "(1) 

Once  the  component  models  are  selected,  drive  train  models  can  be 
synthesized  and  prepared  for  simulation.  No  discussion  will  be  made 
concerning  the  relative  merits  of  digital  simulation  compared  to  analog 
or  analog-digital  (hybrid)  simulation.  The  simulation  literature  is 
replete  with  these  discussions.  Digital  simulation  was  decided  at  the 
outset. 

A  major  part  of  this  study  has  been  concentrated  on  writing  computer 
programs  for  prime  movers,  propellers,  and  transient  operations.  These 
programs  were  developed  to  investigate  various  control  schemes  for  marine 
power  plants. 

The  gas  turbines  selected  for  simulation  are  the  Pratt  and  Whitney 
FTS+A-2  and  the  General  Electric  LM2500-A.  A  medium  speed  and  a  low  speed 
diesel  were  also  selected.  The  medium  speed  diesel  selected  is  the 
Mirrless  KV  Major  12;  the  low  speed  diesel  selected  is  the  Burraeister 
and  Wain  7K98FF.  The  main  reasons  for  selecting  these  particular  engines 
were  the  availability  of  suitable  data  at  reasonable  power  levels. 


1.  GAS  TURBINE  MODELS 

1.1  Introduction 

Gas  turbines,  as  applied  to  ship  propulsion,  usually  connote 
the  inclusion  of  a  power  turbine  or  "free"  turbine,  which 
provides  the  driving  torque.  The  power  turbine  is  driven  by 
the  hot  exhaust  gases  of  the  "gas  generator."  Thus  when 
marine  engineers  refer  to  gas  turbines,  they  usually  include 
the  combination  of  gas  generator  and  power  turbine  in  their 
meaning. 

There  are,  of  course,  exceptions  to  the  statement  con- 
cerning separate  power  turbines.  The  Koelin  Class  escort 
frigates  built  for  the  German  Federal  Navy  are  driven  by 
gas  turbines  which  do  not  have  a  separate  power  turbine. 
The  gas  generator  is  coupled  through  a  clutch  to  the  drive 
train. (2) 

1.2  Aircraft  Gas  Turbine  Models 

During  the  literature  survey  part  of  this  study,  several 
gas  turbine  models  useful  for  aircraft  propulsion  simulation 
were  found.  In  Reference  (3),  for  example,  Cottington  pre- 
sents a  simulation  study  of  a  tx^o-spool  turbojet  engine,  and 
also  covers  the  design  of  a  digital  three-terra  controller 
using  digital  simulation  methods.  In  trying  to  determine  an 
"optimum"  set  of  gains,  he  treated  the  high  pressure  rotor 
system  as  a  critically  damped  second  order  system.  Values 
of  these  optimized  gain  settings  for  varying  values  of  high 
pressure  rotor  speeds  are  presented. 

Mueller  has  developed  a  single  spool  gas  turbine  model 
from  one-dimensional  flow  relationships. (*0  He  assumes  that 
the  only  important  dynamics  are  those  of  the  rotating  compressor- 


1.  GAS  TURBINE  MODELS  (Cont'd) 

1*2  Aircraft  Gas  Turbine  Models  (Cont'd) 

turbine  spool.  Mueller  justifies  this  by  noting  that  gas 
disturbances  travel  through  the  engine  at  the  speed  of  sound, 
whereas  a  spool  acceleration  or  deceleration  takes  some  seconds. 
No  results  of  his  simulations  are  presented. 

In  References  (5)  and  (6),  models  are  derived  from  a 
linear  analysis  of  gas  turbine  dynamics.  The  model  presented 
in  the  latter  reference  is  highly  detailed  and  is  developed 
using  non-dimensionalized  parameters.  To  use  this  model, 
one  must  have  access  to  the  machinery  performance  maps.  The 
prominent  feature  of  this  model  is  that  it  assumes  stepwise 
linearity.  The  compressor  performance  map  is  utilized  to 
obtain  the  coefficients  in  the  linearized  equation  at  each 
operating  point. 

In  References  (7)  and  (8),  Saravanamuttoo  and  Fawke  dis- 
cuss the  application  of  both  analog  and  digital  computer  methods 
of  simulating  gas  turbine  dynamic  performance.  They  approached 
the  problem  from  the  viewpoint  of  engineering  thermodynamics, 
using  the  normal  compressor  and  turbine  characteristics. 
Figure  (1-1)  shows  the  results  obtained  by  Saravanamuttoo  and 
Fawke  for  a  single-spool  turbojet  simulation,  both  digital 
and  analog,  along  with  the  results  of  actual  engine  tests. 
It  can  be  seen  that  the  acceleration  predicted  by  both  simu- 
lations was  much  faster  than  the  actual  engine.  An  examination 
of  the  test  results  on  the  engine  showed  that  there  was  a  lag 
in  the  rise  of  turbine  entry  temperature  following  a  fuel  change. 
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1.  GAS  TURBINE  MODELS  (Cont'd) 

1.2  Aircraft  Gas  Turbine  Models  (C.>r.t*d) 

Saravanamuttoo  and  Fawke's  model,  however,  had  assumed 
Immediate  response,  although  they  realized  there  would  be 
some  lag  due  to  the  dynamics  of  the  fuel  injection  system. 
They  decided  to  simulate  this  lag  by  postulating  that  after 
a  sudden  change  in  fuel  flow,  the  combustion  efficiency  drops 
instantaneously  and  then  recovers  exponentially.  This  drop 
in  efficiency  was  made  proportional  to  the  percentage  increase 
in  fuel  flow.  The  constant  of  proportionality  was  determined 
such  that  the  turbine  entry  temperature  before  and  after  the 
transient  was  unchanged  from  the  steady  state  values*  The 
tirce  constant  for  tha  recovery  of  combustion  efficiency  was 
taken  to  be  inversely  proportional  to  the  air  flow  through 
the  combustor,  i.e., 

2T=  21  .  &s 

*£  -  time  constant  for  recovery  of  combustion  efficiency 

£0  -  time  constant  at  design  speed 

W0-  air  flow  at  design  speed 

W  -  air  flow 

To  achieve  flow  compatability  in  their  digital  model, 
Saravanamuttoo  and  Fawke  added  lags  in  the  propagation  of 
pressure  signals  through  the  engine.  The  lags  were  intro- 
duced by  considering  the  engine  components  to  be  separated 
by  specified  volumes;  these  volumes  produced  the  desired 
capacitive  effect. 


1.  GAS  TURBINE  MODELS  (Cont'd) 

1»3  Marine  Gas  Turbine  Models 

During  the  early  parts  of  the  literature  survey,  a  Russian 
title  on  analog  and  digital  models  of  marine  gas  turbines 
was  found  in  a  gas  turbine  newsletter.  This  book,  listed 
as  Reference  (9),  was  written  by  Bichayev  and  published  by 
the  Shipbuilding  Press  of  Leningrad  in  1969.  Although  it 
is  still  untranslated,  the  book  appears  to  be  a  major  work 
in  the  field  of  marine  gas  turbine  simulation,  judging  from 
the  many  analog  computer  diagrams,  compressor  maps,  and 
equations.  The  cost  of  having  it  translated  (over  $900) 
was  prohibitive. 

An  unedited  rough  translation  of  a  Russian  work  by  Lurye 
discusses  mathematical  modeling  of  a  marine  gas  turbine  engine 
with  a  real  fuel  apparatus. (10)  Lurye's  model  is  presented 
in  the  form  of  linearized  differential  equations.  The  equations 
describing  Lurye' s  model  are  reproduced  in  Appendix  A. 

Using  Saravanamuttoo » s  technique,  one  can  extend  his 
model  to  a  gas  turbine  having  a  free  turbine.  This  is  also 
presented  in  Appendix  A. 
1.^  Torque-Speed.  Characteristics  Model 

C.  J.  Rubis  has  developed  a  method  based  on  the  free 
turbine  torque-speed  characteristics  to  model  a  gas  turbine's 
action. (11, 12)  The  free  turbine  torque-speed  characteristics 
can  be  derived  from  the  steady-state  power-speed  characteristics 
which  show  lines  of  constant  SFC.  The  corresponding  curves 
for  the  Pratt  and  Whitney  FT^A-2  and  General  Electric  LH2500-A 


1.  GAS  TURBINE  HOTELS  (Cont'd) 

1A  Torque-Speed  Characteristics  Model  (Cont'd) 

are  given  in  Appendix  A.  The  torque-speed  characteristics 
have  been  used  here  as  the  model  selected  for  simulating  the 
performance  of  the  marine  gas  turbine. 

The  unrefined  model,  used  to  describe  the  transient 
performance  of  the  turbine  operates  in  the  following  manner. 
Referring  to  Figure  (1-2),  assume  that  the  initial  steady-state 
operating  point  is  at  A  operating  at  torque  TQ  and  speed  N^, 
and  assume  that  the  engine  is  to  accelerate  to  point  B.  In 
a  time  interval  At,   fuel  is  increased  by  A  F.  This  increase 
in  fuel  produces  an  instantaneous  change  in  torque  A  T  such 
that  the  new  torque  is  TQ  +  Z^T.  This  increase  in  torque  will 
result  in  a  change  in  power  turbine  rotation  A  N  along  the 
constant  fuel  rate  line  FQ  +4F.  The  new  starting  point  for 
the  next  time  interval  is  at  the  intersection  of  the  F0+AF 
constant  fuel  rate  line  and  the  new  rotational  speed  L+ANi 
The  dotted  lines  illustrate  the  paths  taken  during  the  time 
steps,  while  the  heavy  line  is  the  torque  path. 

1.5  An  Improved  Model 

The  methodology  of  this  model  is  not  new.  The  torque- 
speed  characteristics  have  been  used  in  the  past  to  model 
diesel  engines.  The  idea  of  obtaining  instantaneous  torque 
for  a  change  in  fuel  rate  is  much  closer  to  the  real  thing  in 
a  diesel  than  it  is  in  a  gas  turbine.  Comparisons  with  engine 
tests  have  shown  that  the  predicted  values  of  torques  obtained 
from  the  gas  turbine  characteristics  are  higher  than  the 
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1.     GAS  TtTRBItlE  MODELS  (Cont'd) 

1.5  An  Improved  Model  (Cont'd) 

measured  values*  The  differences  in  torque  and  the  resulting 
differences  in  speed  were  measured  by  Rubis  who  conducted 
engine  tests  with  the  FT^A-2.  He  also  used  the  torque-speed 
curves  in  a  simulation  between  the  same  two  power  levels. 
Figure  (1-3)  shows  how  the  simulation  results  in  a  higher 
acceleration  because  of  the  higher  torques.  This  graph  should 
be  compared  with  figure  (i-i). 

It  is  evident  from  the  comparison  of  figures  (i-J)   and  (>-'$ 
that  what  is  needed  to  improve  this  model  is  a  correction 
factor,  which  could  have  a  physical  explanation  similar  to 
that  offered  by  Saravananuttoo  and  Fawke  for  their  model; 
namely  a  reduction  in  combustion  efficiency.  The  results  of 
other  transient  tests  conducted  with  the  FT^A-2  were  presented 
by  Rubis  in  Reference  (13).  In  this  paper,  Rubis  presents 
correction  factors  to  the  torque-speed  curve  model  for  engine 
accelerations.  This  curve  is  reproduced  herein  as  figure  (1-^). 
One  can  observe  that  the  curve  for  accelerations  from  idle  can 
be  expressed  as  a  decaying  exponential  similar  to  the  proposal 
by  Saravanarauttoo  and  Fawke,  viz., 

Ke  =  .3  +  .7  U-e'*7"5*  ) 
This  expression,  which  creates  a  form  of  delay  or  lag  in  the 
transient  torque  produced  by  the  model,  was  included  in  the 
simulations  in  this  study.  It  should  be  clear  that  this  lag 
is  applied  to  the  transient  torque  obtained  from  the  model 
and  not  to  the  whole  torque. 
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K     =  TRANSIENT   ENGINE  TORQUE 
E  ="   STEADY  STATE    ENGINE  TORQUE 


4  6  8 

TIME  .  sec 
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CORRECTION  FACTORS  TO  STEADY  STATE  ENGINE  TORQUE 
MODEL 


Fig.  (l-<0 
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!•  GAS  TURBINE  MODELS  (Cont'd) 

1.5  An  Improved  Model  (Cont'd) 

The  model  appears  to  be  acceptable  for  ship  propulsion 
studies  because  it  is  simple;  it  also  permits  loading  other 
characteristics  of  the  ship  drive  train  into  a  medium  sized 
computer  without  exceeding  its  core  size. 
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2.     DIESEL  ENGINE 

2.1  Introduction 

The  traditional  use  of  steam  power  plants  to  propel 
large  American  built  ships  has  created  somewhat  of  a  void 
in  this  country *s  technical  literature  on  diesel  ship  pro- 
pulsion. The  successful  use  of  diesel  engines  for  ship 
propulsion  by  foreign  shipbuilders,  on  the  other  hand,  has 
been  an  established  fact  long  before  the  first  digital  com- 
puter was  ever  built,  and  long  before  analog  and  digital 
simulations  became  established  as  accurate  techniques  for 
modeling  and  testing  ship  propulsion  plants  prior  to  their 
construction.  Consequently  there  has  been  very  little  demand 
from  industry,  either  foreign  or  domestic,  for  diesel  plant 
simulations.  However,  some  published  data  were  obtained 
and  used  in  this  study.  Since  the  use  of  high  speed  diesels 
has  been  fairly  restricted  to  small  craft,  the  simulations 
have  been  limited  to  medium  and  low  speed  diesels. 

2.2  Presentation  of  Models 

Woodward  presents  in  Reference  (1^)  a  simple,  straight- 
forward treatment  of  the  diesel  engine  as  applied  to  ship 
propulsion.  A  widely  used  expression  given  in  Reference  (1*0 
and  other  engine  textbooks,  is  the  formula  for  indicated 
horsepower  (IHP): 

ffiP  =  Ho  (HP> 

P  =  indicated  mean  effective  pressure  (lbs/in  ) 

L  =  stroke  (ft) 

A  =  total  piston  area  (in  )  for  a  multi-cylinder  engine 

RPM 

N  =  RPM  (two-stroke  engine),  •—  (four-stroke  engine) 
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2.  DIESEL  ENGINE  (Cont'd) 

2.2  Presentation  of  Models  (Cont'd) 

This  power  is  that  produced  in  the  engine  cylinders,  and  is 
greater  than  that  available  at  the  output  shaft  because  of 
friction  and  other  internal  loads  (e.g.,  water  pumps,  fuel  pumps, 
turbochangers,  etc.).  What  is  available  at  the  output  shaft 
is  termed  brakehorsepower  (BHP).  The  ratio  of  brakehorsepower 
(BHP)  to  indicated  horsepower  (IKP)  is  the  mechanical  efficiency 
(77  )  of  the  engine,  i.e., 

»     =  £2£ 

"m    MP 

or  BHP  =  V   PUN 
3305o~ 

The  product  N  P  is  called  the  brake  mean  effective  pressure 
m 

(BMSP).  The  indicated  mean  effective  pressure,  P,  can  be 

directly  measured  and  is  proportional  to  the  energy  released 

by  combustion  in  the  cylinder.  This  energy  is  proportional 

to  the  amount  of  fuel  injected  in  each  power  stroke.  Most 

diesels  utilize  positive  displacement  pumps  for  fuel  injection. 

The  amount  injected  per  power  stroke  for  this  type  of  pumping 

system  should  be  the  same.  The  stroke,  L,  and  the  total  piston 

area,  A,  are,  of  course,  fixed  for  a  particular  diesel.  If 

one  assumes  constant  mechanical  efficiency,  then  it  can  be 

seen  that  for  a  particular  fuel  setting,  BHP.  is  a  linear 

function  of  N,  i.e., 

BHP  =  KN 

K  =  N  PXA    for  a  particular  fuel  injection  setting 
33000 
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2.  DIESEL  ENGINE  (Cont'd) 

2.2  Presentation  of  Models  (Cont'd) 

Since  rotational  power  is  the  product  of  torque,  £,  and  speed, 
Co,   it  then  follows  that  for  a  two  stroke  engine, 
co  =  2  7T  N   (rad/rain) 
Zco  =    %  2  ff  n   (ft-lb/min) 
2^=  K_  K  as  previously  defined 

i.e.,  torque  is  constant  for  particular  fuel  injection  setting. 
The  above  is  known  as  the  ideal  torque-rpm  and  power-rpm 
diesel  relationships,  and  are  graphically  illustrated  in 
Figure  2-1  and  2-2. 

The  important  assumption  of  this  model  is  the  one  concern- 
ing constant  mechanical  efficiency.  How  close  does  this 
approach  the  actual  case?  In  Reference  (15) >  Pounder  devotes 
an  entire  chapter  to  the  Burraeister  and  Wain  Engines  and  states 
a  mechanical  efficiency  of  91$  in  his  discussion  on  B  &W  low 
speed  diesels.  (Note:  this  figure  was  confirmed  by  a  representa- 
tive  of  the  Burmeister  &  Wain  Corporation  on  a  recent  visit  to 
MIT). 

It  was  indicated  before  that  the  indicated  mean  effective 
pressure  is  proportional  to  the  fuel  injected  rate.  An  ex- 
tension of  the  ideal  model  then  is  to  say  that  torque  is  not 
only  constant  for  a  given  fuel  setting,  but  also  proportional 
to  that  fuel  setting.  In  other  words,  torque  is  a  linear 
function  of  fuel  setting.  Obviously,  other  extensions  can 
be  made  such  as  relating  torque  to  fuel  rack  position  or  to 
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2.  DIESEL  ENGINE  (Cont'd) 

2.2  Presentation  of  Models  (Cont'd) 

the  movement  of  some  fuel  control  actuator, 

Lewis,  Lecourt,  and  Scoville  used  a  variation  of  the 
ideal  diesel  engine  in  Reference  (16).  The  torque  developed 
within  the  engine  was  designated  Q   and  the  friction  torque 
Q»  •  The  output  torque  is  the  difference  between  the  two 
torques.  Qp„  was  postulated  as  linearly  related  to  fuel 
rack  position,  z,  i.e., 

and  the  friction  torque  was  taken  as  linearly  related  to  engine 
speed, 

Qfe  =Ab+ArN. 
It  can  be  seen  that  this  model  does  not  assume  constant 

mechanical  efficiency. 

In  Reference  (17)  >  Garvey  presents  a  linearised  differen- 
tial equation  describing  the  motion  of  a  diesel  engine.  He 
assumed  the  output  torque  developed  by  the  engine  to  be  pro- 
portional to  governor  position,  i.e., 
T  =  Cz 

T  -  change  in  torque 
z  -  governor  position 
c  -  governor  gain 
Parker  uses  the  same  model  with  a  delay  (18), 
T  =  ce   z 
T  =  Torque 
c  *  governor  gain 
z  =  governor  position 
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2.  DIESEL  ENGINE  (Cont'd) 

2.2  Presentation  of  Models  (ContM) 

s  =  differential  operator 
f  =  time  delay 

2.3  Low  Speed  Diesel  Model 

The  ideal  diesel  engine  model  in  conjunction  \^th  a  delay 
in  the  form  of  an  exponential  lag  was  chosen  as  the  method 
to  simulate  a  Burmeister  &  Wain  7K98FF  low  speed  diesel. 
In  Reference  (19) »  Keulengracht  presents  the  results  from 
test-bed  trials  of  a  7K98FF.  From  these  results,  and  assuming 
the  ideal  diesel  characteristics,  one  can  get  a  brakehorse- 
power-speed  curve  and  from  this  curve,  derive  the  flat  torque 
speed  characteristics  for  constant  fuel  settings.  These 
curves  are  all  included  in  Appendix  B. 
ZA     Medium  Speed  Diesel  Model 

Reference  (20)  describes  the  development  of  high  per- 
formance medium  speed  diesel,  the  Kirrless  KV  Major  12. 
Since  the  fuel  map  for  this  diesel  is  included  in  this  article, 
it  was  selected  as  the  diesel  to  model*  From  the  fuel  map, 
and  using  the  elementary  relations  presented  earlier,  one 
can  derive  the  engine* s  torque  speed  characteristics  for  con- 
stant fuel  settings.  In  the  process  of  drawing  these  charac- 
teristics, it  was  observed  that  torque  could  be  approximated 
quite  well  as  being  proportional  to  fuel  setting.  Thus 
the  expression, 

T  =  ICx  (torque) 
x  =  fuel  setting 
K  =  proportionality  constant 

was  used  in  the  simulations. 
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2.     DIESEL  EIIGIHE  (Cont'd) 
2.5     Delays 

As  in  the  gas  turbine,  the  diesel  engine  does  not 
instantaneously  produce  a  change  in  torque  with  an  increase 
in  fuel  rate.  For  a  four-stroke  cycle  multi-cylinder  diesel 


engine,  Parker  recommends  a  value  of  rf-  when  N  equals  engine 

N 


11 
N 

rpm;  this  is  an  approximation  for  diesel  engines  with  two 
or  more  cylinders  firing  per  revolution  . (21)  For  the  more 
general  case,  Reference  (16)  recommends 

"Re  NRe      NRep 

Np  :  designed  engine  speed 

p  :  number  of  cylinders  firing  per  revolution 
2.6  Discussion 

The  torque-fuel  setting  characteristics  were  found  to 
be  linear  for  the  medium  speed -diesel  but  not  quite  linear 
for  the  low  speed  engine.  The  characteristics  for  the  low- 
speed  diesel  were  therefore  stored  in  the  computer.  For  each 
diesel,  a  change  in  fuel  setting  results  in  a  change  in  torque. 
This  torque  was  tempered  by  a  delay  in  the  form  at  a  decaying 
exponential.  This  method  is  much  simpler  than  that  employed 
for  the  gas  turbine  models,  since  diesel  torque  is  independent 
of  speed  in  the  engine fs  operating  range. 

For  ship  propulsion  studies  being  considered  here,  it  was 
not  considered  necessary  to  seek  a  model  which  would  produce 
the  oscillatory  torque  that  is  characteristic  of  internal  com- 
bustion engines  (cf.  Reference  16). 


20 


3.  CLUTCHES  AND  HYDRQDYNAKEC  TORQUE  CONVERTERS 
3.1  Introduction 

The  selection  of  a  fixed  pitch  propeller  and  a  non- 
reversing  prime  mover  will  necessarily  mean  the  incorporation 
of  a  clutch  or  torque  converter  in  the  drive  train.  Gas 
turbines,  for  example,  are  limited  to  undirectional  movement; 
and  some  diesels,  on  the  other  hand,  are  designed  for  operation 
in  two  directions.  It  v&ll  be  observed  later  that  clutch 
engagement  and  disengagement  times  play  a  critical  part  in 
the  emergency  stopping  of  a  ship.  The  selection  of  a  proper 
clutch  for  a  drive  train  thus  is  no  trivial  matter.  In  the 
same  vein,  the  selection  of  a  suitable  clutch  model  to  simu- 
late the  emergency  stopping  of  a  ship  is  important  in  order 
to  avoid  erroneous  results. 
3«2  FaJk  Marine  Airflex  Clutch 

Figure  3-1  shows  the  principal  part  of  the  Airflex  clutch. 
The  clutch  is  pneumatically  operated  drum-type  dry  friction 
clutch  consisting  of  an  annular  rubber/fabric  air  gland  bonded 
to  an  outer,  flanged  steel  rim  and  lined  with  friction  shoes 
which  engage  an  inner  iron  drum  when  the  gland  is  inflated. 
A  newer  air-cooled  version  utilizes  extruded  aluminum  shoe 
blocks  with  axial  air  passages,  to  which  the  friction  brake 
block  liners  are  riveted  to  insulate  the  rubber  gland  from 
the  heat  generated  at  the  friction  surface  (cf.  Reference  22). 

Rubis  was  able  to  obtain  from  the  Falk  Corporation  the 
torque  expressions  for  the  clutches. (11)  The  torque  capacity 
was  found  to  be  directly  related  to  the  number  of  brake  liners 
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FALK  MARINE  AIHFLEX  CLUTCH 


Fig.  (3-D 
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3*  CLUTCHES  AND  HYDRODYNAMIC  TORQUE  CONVERTERS  (Cont'd) 
3*2  Falk  Marine  Airflex  Clutch  (Cont'd) 

and  tho  net  air  pressure  used  to- engage  the  clutch.  For  a 
35-inch-diameter  clutch  with  eight  liners,  the  maximum  clutch 
torque  capability,  Q  ,  is  given  by  the  expression 

%    =2275Pnet   (lb-ft> 

Pne+  =  ne"t  sir  pressure  activating  the  clutch  (psi) 

One  of  the  important  features  to  consider  in  this  clutch 

is  that  the  inward  expansion  of  the  clutch  glands  during 

engagement  is  opposed  by  the  centrifugal  force  of  the  rotating 

clutch  assembly.  The  retarding  force  has  been  transformed 

to  an  equivalent  air  pressure,  pc,  and  is  expressed  by  the 

relation 

pc  =  5  +  (5.9  x  10"5)Ng2     (psi) 

Ng  =  speed  of  the  frame  housing  the  clutch  glands  (rpm) 
Rubis  states  that  typical  inflation  and  deflation  rates  are 
5  psi  per  second  and  30  psi  per  second  respectively;  and  a 
typical  supply  pressure  is  150  psi. 
3»3  SQ500  Synchroclutch  (Dental  Clutch  Type) 

The  SQ500  Synchroclutch  is  manufactured  by  the  Philadelphia 
Gear  Corporation.  Literature  on  the  description  and  operation 
was  furnished  by  the  manufacturer  and  is  included  in  Appendix  C. 
This  particular  clutch  is  recommended  for  gas  turbines  in  the 
25000  HP  at  36OO  rpm  range  (e.g.  FT^A-2  &  LK2500).   For  a 
typical  application,  the  torque  history  for  the  synchronizing 
(friction)  clutch  may  be  approximated  by: 
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CLUTCHES  AMD  HYDRODYNAHIC  TORQUE  CONVERTERS  (Cont'd) 
3.3  SQ500  Synchroclutch  (Dantal  Clutch  Type)  (Cont'd) 

T  =  T*(Wt/26> 

T  =  torque  at  tine  t  (ft-lb) 

t  =  time  (sec) 

T1  =  11,600  (ft-lb) 
Both  IT  and  the  time  constant  may  be  changed  to  suit  specific 
operating  conditions. (23)  Other  information  and  details  on 
this  clutch  are  included  in  Appendix  C. 

Another  clutch  which  operates  in  a  similar  manner  to  the 
SQ500  is  the  BIH-I}ynetic  synchronizing  clutch  manufactured 
by  Koppers  Co,  Inc.  Although  no  torque-time  history  expression 
was  provided  by  the  manufacturer,  information  was  provided 
on  the  operation  and  design  of  the  clutch.  This  is  included 
in  Appendix  C  for  information  purposes  only. 
3.^  Hydraulic  Transmissions 

One  of  the  problems  inherent  in  a  plant  driven  by  a  diesel 
or  gas  turbine  prime  mover  with  a  fixed  pitch  propeller  that 
is  not  found  in  a  comparable  steam  plant  is  that  the  former 
prime  movers  have  minimum  speeds  at  which  they  can  be  operated. 
The  Coast  Guard  Ship  Hamilton  (WPG-715)*  for  example,  is  powered 
by  two  FT^-A-6  gas  turbines  and  two  Fairbanks  Morse  Model 
38TD-1/8  turbo-supercharged  opposed-piston  marine  diesel  engines 
driving  two  controllable  pitch  propellers.  With  the  propeller 
at  full  pitch,  and  the  diesel  at  idle,  the  vessel  will  attain 
a  speed  of  approximately  12  knots,  and  18  knots  with  full  pitch 


3.  CLUTCHES  AND  HYDRODTiTAHIC  TORQUE  COIWERTKRS  (Cont'd) 
3»^  Hydraulic  Transmissions  (Cont'd) 

and  turbine  at  idle.  (2*1 )  Richardson  describes  a  method 
whereby  the  Falk  clutch  can  be  partially  inflated  thereby 
permitting  continuous  operation  or  maneuvering  at  shaft  speed 
below  that  corresponding  to  minimum  turbine  (or  diesel)  idle. (25) 
For  some  applications,  however,  the  heat  rate  may  be  prohibitive 
for  continually  doing  this,  and  its  effect  on  the  life  of  the 
clutch  is  not  clear. 

In  Reference  (26),  Rowen  describes  a  hydraulic  trans- 
mission system  that  permits  propeller  speeds  in  the  low  rpms 
with  a  diesel  or  gas  turbine.  The  system  has  a  hydraulic 
reversing  gear  which  incorporates  a  fluid  coupling  for  ahead 
rotation  of  the  propeller  and  a  reversing  torque  converter  for 
astern  operation.  In  addition,  a  positive  drive  clutch  is 
included  for  normal  operation  in  the  ahead  direction  at  sea. 
Figure  3-2  is  a  cross  section  of  the  hydraulic  transmission 
system.  The  torque  characteristics  on  a  "per  unit"  basis  of 
the  two  hydraulic  elements  are  presented  in  figure  3-3 •* 
Speed  ratio,  as  defined  by  Rowen,  is  the  ratio  of  the  output 
speed  to  the  input  speed  (positive  for  ahead  rotation,  negative 
for  astern  rotation).  He  assumed  that  the  hydraulic  element 
efficiencies  were  constant  over  the  speed  ratio  range,  and 
defined  the  output  torque  by  the  relation, 

T   =Q»x|!2  xF 
c    *of    pm 

*  "Per  Unit"  is  defined  as  the  actual  value  divided  by  the  design  point  value. 
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Fig.  (3-3) 
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3.  CLUTCHES  AND  HYDRODYNAHIC  TORQUE  CONVERTERS  (Cont'd) 
3*^  Hydraulic  Transmissions  (Cont'd) 

T   =  converter  or  coupling  output  torque 
Q  »  =  oil  flox*  to  converter  or  coupling 
N   =  input  speed  of  prime  mover 
F   -a  function  of  the  speed  ratio  (Fig.  3-3) 
Figure  3-4  is  a  block  diagram  of  the  hydraulic  transmission 
simulation  used  by  Rowen.  No  numbers  were  given  by  him  for 
the  various  parameters  involved  in  the  simulation, 
3*5  Discussion 

The  Falk  clutch  and  the  Philadelphia  Gear  synchroclutch 
were  selected  for  inclusion  in  the  simulation  because  of  their 
simplicity  and  because  of  their  wide  spread  use  in  ship  pro- 
pulsion trains.  Since  no  expression  for  clutch  disengagement 
was  provided  by  Philadelphia  Gear,  it  was  assumed  that  the 
clutch  torque  decreased  as  an  exponential  decay  with  a  time 
constant  of  5  seconds. 
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PROPELLER  MODELS 

^.1  Introduction 

The  propeller  simulation  proved  to  be  the  most  challenging 
aspect  of  this  study.  An  excellent  study  of  propellers  under 
conditions  of  reversing  was  written  by  Miniovich  and  is  listed 
as  Reference  (28).  Miniovich  questions  the  validity  of  using 
propeller  characteristics  measured  under  equilibrium  conditions 
in  the  simulation  of  backing  and  reversing  of  ships.  To  answer 
this,  he  performed  reversing  experiments  with  models  in  a  towing 
tank  and  measured  the  propeller  characteristics.  These  char- 
acteristics were  found  to  be  in  agreement  with  those  measured 
in  equilibrium-  or  quasi-stationary  tests.  The  only  notable 
discrepancy  appeared  to  lie  in  the  zero  thrust  region.  From 
his  equations  of  motion,  Miniovich  also  calculated  the  values 
of  basic  variables  that  characterize  the  reversing  process 
(e.g.,  time  required  to  stop  the  ship,  time  required  to  stop 
the  propeller,  and  maximum  values  of  thrust  and  torque)  and 
found  these  to  differ  from  experimental  data  by  a  maximum 
of  10  percent.  It  thus  appears  that  one  can  use  the  steady 
state  propeller  characteristics  in  a  simulation  and  get  a 
reasonable  accuracy. 

k,Z    Propeller  Characteristics 

The  traditional  definition  of  the  propeller  advance 
coefficient  is  the  familiar  expression, 

J  =  Ik 

ND 

V.  =  propeller  speed  of  advance  (ft/sec) 

A 
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^.  PROPELLER  MODELS  (Cont'd) 

*K2  Propeller  Characteristics  (Cont'd) 

N  =  propeller  rate  of  rotation  (rev/ sec) 

D  =  propeller  diameter  (ft) 

This  coefficient  can  be  derived  from  a  dimensional  analysis 

for  propellers  (cf.  Reference  29)  •  It  is  obvious  that  relating 

fixed  pitch  propeller  coefficients  to  J  might  prove  to  be 

unwieldy  for  ship  propulsion  studies  since  J  becomes  unbounded 

as  the  shaft  goes  from  positive  rotation  (providing  thrust 

to  move  the  ship  ahead)  to  negative  rotation  (providing  thrust 

to  move  the  ship  astern),  Miniovich  solved  this  dilemma  by 

defining  the  propeller  characteristics  in  terms  of  both  J 

and  its  inverse,  viz.,  — -   _,      ,.  .   .      ,   ^ 

V&.  The  result  is  two  graphs  of 

propeller  characteristics  instead  of  one;  when  either  one 

of  the  coefficients  tends  to  large  values  he  changes  to  the 

alternate  definition. 

At  the  Ninth  International  Towing  Tank  Conference  (1957) » 

Stephanson  demonstrated  how  the  J  and  the  l/j  propeller 

characteristics  could  be  combined  with  a  new  term,  called 

the  modified  advance  coefficient  defined  by, 


Cv  ~ 

A 
where  all  the  terms  are  defined  as  before. (30)  (Miniovich 

ascribes  a  similar  expression  to  Kalmakov  in  19^0).  It  can 

be  seen  that  J'  will  not  become  unbounded  as  N  or  V,  approach 

zero  values. 


31 


**•  PRQPELIER  MODELS  (Cont'd) 

^.2  Propeller  Characteristics  (Cont'd) 

In  Reference  (31)  Baker  and  Patterson  discuss  Miniovich's 
work  and  present  his  propeller  characteristics  using  the  modi- 
fied advance  coefficient.  Other  forms  of  advance  coefficients 
are  offered  by  the  two  authors  to  make  the  task  of  simulating 
a  propeller  easier  than  as  performed  by  Miniovich. 
^•3  Modeling  the  Wageningen  B-Screw  Series 

In  1969 t  Van  Lammeren,  Van  Manen,  and  Oosterveld  presented 
a  new  technique  of  reproducing  the  propeller  characteristics 
of  the  Wageningen  B-Screw  Series. (32)  Th9  nondimensionalized 
thrust  and  torque  coefficients  are  presented  in  the  form  of 
a  Fourier  series.  The  series  is  given  as  a  function  of  the 
hydrodynamic  pitch  angle,  A9   defined  by, 

A    =  arctan  (  _A_   ) 

where  V.jN*  and  D  are  as  previously  defined.  The  nondimen- 

sionalized  thrust  coefficient,  C„*,  and  torque  coefficient, 

C  *,  are  then  given  by 

C  *  =  M_ 
1     ^   A(K)  cos  £K  +  B(K)  sin  AK 

10  C  *  «  M 

q    ^~  A(K)  C03  y8K  +  B(K)  sin  f&L 

Once  the  coefficients  are  obtained,  thrust,  T,  in  pounds, 
and  torque,  Q,  in  ft-lbs  can  be  calculated  using, 

T  =  ^*{if   ^A2  +  (0*7ffN  D  /\  4-  D2j 
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fc.     PROPELLER  MODELS  (Cont'd) 

^•3     Modeling  the  Wageningen  B-Screw  Series  (Cont'd) 

Q  =  cq*  f  if[\2     +  (0-7  N  v  D  'fJSz^l 
where 

/>  =  fluid  density 

V.,  N»  D  as  previously  defined 

The  factor  "10"  in  the  expression  for  nondimensionalized 

torque  coefficient  apparently  through  omrnission,  does  not 

appear  in  two  different  published  versions  of  Reference  (32). 

Simulations  were  performed  without  this  factor  and  resulted 

in  excess  propeller  torques  (10  tiroes  that  predicted  from 

hand  calculations).  Figure  (kQ)   of  Reference  (32)  was 

reproduced  by  the  IEM  1130  computer  without  the  factor  of 

10  in  the  expression.  Examination  of  Figure  (^0)  revealed 

that  the  ordinate  is  labeled  C  *  and  -10  C  *,  the  minus  sign 

T         q 

being  introduced  to  avoid  overlapping  of  curves.  It  was  thus 
deduced  that  the  expression  for  C  *  should  include  a  factor  of  10, 

q 

The  values  of  the  coefficients  A(K)  and  B(K)  for  a  series 
of  20  terras,  i.e.  K  =  1-20,  are  presented  in  Reference  (32), 
These  coefficients  are  presented  for  a  wide  range  of  B-Screw 
series.  Figure  ^-1  illustrates  how  closely  a  series  of  20 
terras  matches  the  measured  characteristics  of  a  B  ^-70  screw 
with  a  pitch  ratio  of  1,0,  Figure  ^-2,  furthermore,  illustrates 
how  well  a  series  of  10  terms  approximates  the  measured  char- 
acteristics for  the  same  propeller. 
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A  COMPARISON  OF  A  20   T3KM  FOURIER  SERIES 
•  IT.;  ACTUAL  PROPELLER  CHARACTERISTICS 


Fig.  (4-1) 
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A  COMPARISON  OF  A  10  T,::;.  FOURIER  SERIES 
IsITH  ACTUAL  PROPELLER  CHARACTERISTICS 


Fig.  C*-2) 
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b*     PROPELLER  MODELS  ( Cont  \ d  ) 

^•3  Modeling  the  Wageningen  D-Screw  Series  (Cont'd) 

It  was  probably  noticed  that  the  definition  of  A   does 
not  avert  the  problems  previously  associated  with  the  tradi. 
tional  definition  of  J,  i.e.,  dividing  when  n=o.  It  would 
seem  more  appropriate  to  use 


&    =  arcsine   (V^  +  (0.7  ^  IID)*)^ 

or  0>7   HP       k 

fl  =  arccosine  Qi  c    +   (0.7  jntD))2 

instead  of  the  arctangent  definition  of  A      •  The  IBM 

scientific  subroutine  library  on  the  IBM  1130  digital  computer 

that  was  used  did  not  have  subroutines  for  calculating  arccosine 

or  arcsine.  The  arctangent  subroutine  had  to  be  used.  It 

should  be  also  noted  that  the  IBM  subroutine  for  arctangent 

generates  values  only  between  -  tf£  ^°    +  0L        It  is  the 

2        2  . 

writer's  opinion  that  if  a  numerical  overflow  should  occur 
in  the  argument  for  the  arctangent  subroutine,  that  this  sub- 
routine will  yield  either  ^/2  or  -  ™/2  depending  on  the  sign 
of  the  overflow.  No  numerical  overflow  problems  were  experienced 
using  the  arctangent  subroutine. 
bA    Controllable  Reversible  Pitch  Propellers  (CRPP) 

A  simple  model  of  a  CRPP  was  presented  by  Rowen  in 
Reference  (26).  The  "per  unit"  equation  for  thrust  is  given 

by 

2 
Thrust  =  (pitch  x  rpm) 


**.  PROPELLER  MOIELS  (Cont'd) 

bA     Controllable  Reversible  Pitch  Propellers  (CRPP)  (Cont'd) 
and  the  "per  unit"  equation  for  torque  is  given  by 
Torque  =  £0.3  x  (rpm)  j  +  [0.7   +  (rpw)2  x  (pitch yj 
In  Reference  (33)»  Read  presents  a  study  of  gas  turbines 
driving  a  CRPP,  and  includes  the  torque  and  thrust  characteristics 
for  the  propeller.  Apart  from  the  characteristics,  all  that 

is  known  about  the  propeller  is  that  it  has  '+  blades  and  has 

2 
a  polar  moment  of  inertia  of  7050  lb-ft-sec  .  Wendel  and 

Dunne  apparently  found  a  small  error  in  the  va?.ues  of  modified 
torque  coefficient  given  by  Read  and  applied  a  correction  to 
the  torque  characteristic. (3^)  The  modified  advance  coeffi- 
cient, J',  is  as  previously  defined,  viz., 
J*  =  VA 


(VA2  +  N2D2)f 

The  modified  thrust  coefficient  is  defined  by 

T 
KT   f>  D^V^   +  NV) 

and  the  modified  torque  coefficient  is  defined  by 

KQ'  = 9 

P   D3^2  +  N2D2) 
where  T  =  thrust  in  lbs 
Q  =  torque  ft-lbs 

0  ,  V. ,  N,  D  are  as  previously  defined 
(Note  the  absence  of  7?  in  the  denominator  in  the  expressions 
for  K  '  and  K  *;  $    is  included  in  the  similar  expressions 
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*K  PROPELLER  MODELS  (Cont'd) 

bA     Controllable  Reversible  Pitch  Propellers  (CRPP)  (Cont'd) 
given  by  Van  Lammeren.  Baker  and  Patterson  discuss  the 
"  if  question"  in  Reference  (31)) • 

Figures  (*K3)  and  (Mf)  are  the  CRPP  characteristics,  and 
include  the  Wendel-Dunne  correction. 
^•5  Pitch  Control  System  Model 

The  only  analytical  description  of  a  CRPP  pitch  control 
system  found  was  that  described  by  Read  in  Reference  (33)* 
Pitch,  R,  and  angular  rotation, 0,  are  related  by  the  expression 
R  =  Tan  6 

0^53 

The  pitch-changing  cylinder  has  a  total  capacity  of  ^0  gallons, 
and  the  blade  angle  varies  very  nearly  linearly  with  ram  stroke, 
such  that 

ffi   =  £i-x  &    (Rad/sec) 
dt     57.3   58 

Fo  =  oil  flow  in  gal/sec 

By  manipulating  the  two  equations,  one  can  arrive  at  the 

expression 

dR   =  t   0.0528  +  0.C109  R2Frt 
Ht  ° 

There  are  two  pumps  supplying  the  oil  flow,  FQ.  One  is 

motor  driven  and  the  other  is  shaft  driven.  The  motor 

driven  pump  has  a  maximum  capacity  of  60  gal/rain  and  the 

shaft  driven  pump  has  a  maximum  capacity  of  100  gal/min. 

The  pumps  are  pressure  controlled  such  that  any  alteration 

in  pitch  demand  causes  the  motor-driven  pump  to  provide  60 
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**•     PROPELLER  MODELS  (Cont'd) 

^•5  Pitch  Control  System  Model  (Cont'd) 

gal/rain  and  the  shaft-driven  pump  to  provide  its  maximum 
output  at  the  shaft  speed  in  question.  The  result  is  that 
during  a  pitch  alteration, 

FQ  =  (1  +  0.^35  N)  (gal/sec) 

N   =  shaft  speed  (rev/ sec) 
4-. 6     Discussion 

'Xhe  fixed  pitch  propeller  model  selected  was  the  Fourier 
representation  of  propeller  characteristics  as  presented  in 
Reference  (32).    Only  the  first  ten  coefficients  were  used 
in  the  representation.  The  resistance  characteristics  of  the 
vessel,  and  the  propeller  diameter  are  two  important  factors 
in  selecting  the  prqser  propeller.  Since  the  characteristics 
for  the  B-'4  series  with  an  expanded  area  ratio  of  .?  were 
readily  available,  it  was  decided  to  use  this  type  of  propeller. 
A  plot  of  K  /j  indicated  a  pitch-diameter  ratio  of  1.0  would 
give  the  maximum  efficiency;  hence  this  propeller  was  chosen. 

Read's  controllable  reversible  pitch  propeller  (CRPP) 
was  chosen.  The  propeller  characteristics  as  modified  in 
Reference  (3*0  were  redrawn  for  constant  values  of  modified 
advanced  coefficient  (J').  The  independent  variable  was 
taken  as  pitch  diameter  (PD)  ratio  and  the  dependent  variables 
were  thrust  and  torque  coefficients.  Figure  (4—3)  shows  the 
thrust  coefficient  vs.  PD  ratio  for  J'  =  -0.2.  This  curve 
is  fairly  typical  of  the  others.  The  linear  approximations 
are  indicated.  Figure  (4-4-)  shows  the  torque  coefficients 
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**•     PROPELLER  MODELS  (Cont'd) 
^.6  Discussion  (Confri) 

vs  PD  ratio  for  J'  =  -0.2  with  the  linear  approximations 
used.  This  curve  is  also  fairly  typical  of  the  other 
curves  of  torque  coefficients  for  constant  J*  and  varying 
PD  ratio. 

Read*s  pumping  system  as  described  earlier  was  also 
used  in  conjunction  with  the  CRPP. 
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5.  SHIP  CHARACTERISTICS 

5.1  Wake  Fraction  and  Thrust  Deduction  Factor 

Wake  fraction,  w,  determines  the  velocity  the  propeller 
is  advancing  at  in  relation  to  the  ship  velocity, 

VA  =  (1-w)  V 

V.  =  propeller  speed  of  advance 

V  =  ship  speed 

w  =  wake  fraction 

Thrust  deduction  factor,  t,  indicates  how  much  propeller 
thrust  is  reduced  due  to  hull  interaction, 

Tl  -  To«  (") 

T   =  open  water  propeller  thrust 
ow 

T   =  actual  thrust 

t   =  thrust  deduction  factor 
(For  a  physical  explanation  of  "t"  and  "w",  and  other  useful 
relationships  involving  these  two  terras,  see  Reference  (29) » 
Chapter  VII,  Section  13), 

The  implication  in  the  two  above  definitions  is  that 
they  are  steady  state  quantities.  Thus  at  equilibrium,  the 
actual  thrust  equals  ship  resistance.  Since  t«l.,  the 
definition  shows  how  much  more  thrust  is  required  over  resis- 
tance to  keep  ship  speed  constant. 

Wake  fraction  and  thrust  deduction  factor  are  not,  as 
one  would  expect,  constant  over  a  range  of  speed.  The  variation 
in  t  and  w  depend  on  the  degree  of  interaction  between  propeller 
and  hull.  Harvald  performed  several  towing  tank  tests  to 
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5.  SHIP  CHARACTERISTICS  (Cont'd) 

5»1  Wake  Fraction  and  Thrust  Deduction  Factor  (Cont'd) 

measure  "t"  and  "w"  under  transient  conditions  for  two 
different  hull  shapes. (35)  He  performed  his  towing  tank 
tests  under  conditions  of  constant  model  ship  velocity  and 
varying  shaft  rpm  and  vice  versa.  Although  his  work  pro- 
vided good  insight  into  the  behavior  of  "t»  and  »*w"  under  these 
two  conditions,  it  was  difficult  to  assess  how  his  results 
could  be  applied  for  conditions  of  varying  shaft  rpm  and  ship 
velocity.  How  did  others  treat  this  phenomena?  In  the  solu- 
tion of  his  backing  equations,  Kiniovich  assumed  both  factors 
to  be  equal  to  zero ,(36)  Rubis  assumed  varying  values  of 
"t"  and  "w"  from  two  curves  of  bell  shaped  form. (37)  Read 
considered  wake  fraction  to  be  negligible,  but  formulated 
linear  relations  for  thrust  deduction  factor  for  positive 
and  negative  values  of  propeller  thrust. (38).  Lewis,  Lecount, 
and  Scoville  assumed  thrust  deduction  factor  and  wake  fraction 
to  be  a  function  of  vessel  speed  only;  for  negative  speeds 
thrust  deductions  and  wake  fractions  were  neglected. (39) 
Thus,  it  can  be  seen  that  there  is  no  general  agreement  on 
the  treatment  of  these  two  factors.  The  feeling  appears  to 
be  that  these  two  factors  are  so  small  that  any  errors  in 
their  treatment  will  have  negligible  effects  on  the  total 
simulation.  For  simulations  performed  here,  constant  values 
of  wake  fraction  (w=0.0^)  and  thrust  deduction  factor  (t=0.02) 
were  assumed. 
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5.  SHIP. CHARACTERISTICS  (Cont'd) 

5*2  Added  Mass  of  Ship  and  Propeller 

As  with  thrust  deduction  and  wake  fraction,  the  treat- 
ment of  added  mass  of  ship  and  propeller  is  also  varied. 
The  added  mass  of  the  ship  is  usually  taken  to  be  between 
8»1G$  of  the  displacement.  This  study  used  10$,  Lewis, 
Lecount,  and  Scoville  used  a  different  approach  to  account 
for  added  mass  as  will  be  discussed  later.  The  added  mass 
of  a  fixed  pitch  propeller  is  usually  assumed  to  be  about 
Z%  of  the  inertia.  (cf.,e.g.,  Reference  (^0))  Twenty-five 
percent  was  assumed  here.  There  is  nothing  fixed  about  this 
figure  as  Burrill  and  Robinson  found  the  percentage  addition 
should  be  as  high  as  k6%  for  one  particular  propeller. (^1) 

For  his  CRPP,  Read  proposed  in  Reference  (33)  the  expression, 

I  =  325  +  1050  /R/ 

2 

I  =  added  mass  due  to  entrained  water  (lb-ft-sec') 

R  =  pitch  diameter  ratio 
This  expression  was  used  in  the  simulation  with  the  B  &  W 
7K98FF  low  speed  diesel  -  CRPP  combination. 
5.3  Ship  Resistance 

Accurate  representation  of  ship  resistance  is,  of  course, 
vital  to  the  proper  selection  of  propeller (s).  One  of  the 
basic  assumptions  of  the  popular  "Propeller  Law"  (cf .  References 
(1*0  &  (15)  )  is  that  resistance  is  proportional  to  the  square 
of  velocity  for  ships  where  the  wave  drag  is  small.  For  those 
cases  where  it  is  desired  to  simulate  transients  on  a  parti- 
cular ship,  it  is,  of  course,  more  accurate  to  use  that 
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5<     SHIP  CHARACTERISTICS  (Cont'd) 

50  Ship  Resistance  (Cont'd) 

vessel's  characteristics  as  determined  from  towing  tank  tests. 

Lewis,  Lecount,  and  Scoville  included  acceleration  effects 
in  their  expression  for  ship  resistance, 

n     ss     AD    -rr— 

M   =  mass  of  ship 

V  =  ship  speed 

9 

R_e  «  Rai!S  (V,,)  -   steady-state  resistance 

C.^  =  "added  mass"  coefficient 
AD 

No  values  for  C.n  were  provided  by  the  authors. 

For  simulations  performed  here,  resistance  was  assumed 
to  be  of  the  form, 

R  =  10^.8  V2    (lbs) 

V  -  ship  velocity  (ft/sec) 
5*^  Shaft  Friction 

Read  used  the  expression 

Qf  =  2890  N1*23    (ft-lb) 

N   =  propeller  speed  (assumed  to  be  rev/ sec) 
to  represent  shaft  friction. (33)  Rubis,  on  the  other  hand, 
represented  shaft  friction  characteristics  in  the  form  of  a 
"table- look-up"  in  his  computer  simulations.  These  characteris- 
tics were  approximated  in  simulations  performed  here  by  the 
expressions 

Qf  =  2500  (ft-lb)   0  £h\j  <  25 

Qf     =100/N/     (ft-lb)     /N/  >  25 

N      =  propeller  shaft  rpm 
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6.  RESULTS 


6.1  Introduction 

The  following  maneuvers  were  performed  with  the  drive 
trains  as  indicated: 

FT*JA-2  Gas  Turbine  -  Fixed  pitch  propeller 

1.  Coastdown  from  30  knots 

2.  Accelerating  from  5  knots 

3.  Crashback  from  30  knots 

LJ2500-A  Gas  Turbine  -  controllable  reversible  pitch 
propeller 

1.  Accelerating  from  15  knots 

2.  Crashback  from  30  knots 

Kirrless  K\M-:a.1or  12  -  Fixed  Pitch  propeller 

1,  Crashback  from  9  knots 

Burmeister  &  Wain  7K9&FF  -  Controllable  reversible  pitch 
propeller 

1.  Crashback  from  12  knots 

6.2  Coastdown  from  30  knots  (F'^A°2) 

At  zero  time,  fuel  rate  was  decreased  from  9500  lb/hr  to 
1600  lb/hr  exponentially  with  a  1-second  time  constant. 
Results  of  this  simulation  are  presented  in  figures  (6-1)  - 
(6«<0. 

6.3  Accelerating  from  5  knots  (FT^A-2) 

At  zero  time,  fuel  rate  was  increased  exponentially  with 
a  1-second  time  constant  from  2500  LB/'hr  to  12000  LB/hr. 
Figures  (6-^)   -  (6-11)  are  plots  of  the  results.  A  copy  of 
the  computer  program  used  is  included  in  Appendix  F. 
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6,h     Crashback  from  30  knots  (FF4A-2) 

This  maneuver  was  accomplished  in  three  phases: 

1.  Coastdown*  -  opening  clutch 

2.  engaging  reverse  clutch 

3.  stopping  the  ship 

Results  of  this  simulation  are  plotted  in  figures  (6-12)  - 
(6-17).  Time  for  forward  clutch  disengagement  was  ^-8  seconds. 
The  reverse  clutch  was  engaged  at  11.2  seconds  and  the  ship 
stopping  time  was  ^3»^  seconds.  Head  reach  was  1229  ft. 
A  copy  of  the  computer  program  used  in  phase  3  is  included 
in  Appendix  F» 

6*5    Accelerating  from  15  knots  {IM2500-A) 

A  description  of  the  control  strategy  and  the  controller 
used  is  in  Appendix  G.  Initially  the  ship  is  at  15  knots 
with  a  pitch  ratio  of  1.0  on  the  CRPP.  Pitch  was  increased 
to  1.6  in  8  seconds.  At  that  instant  the  controller  acted 
to  increase  shaft  rpm  to  216.  Results  of  this  simulation 
are  shown  in  figures  (6-18)  -  (6-20 )• 

6»6     Crashback  from  30  knots  (LM2500-A) 

■■■I    ■■!    ■     ■    IP         nil   ■!      ■!■■!     ■     wii      n     »i    iiwi— "*■»      ni-imw  *  i n  mmr         mi  vim 

The  problem  to  consider  here  was  whether  or  not  the  power 
turbine  would  overspeed  while  the  CRPP  goes  through  a  zero 
value  of  torque  coefficient.  This  ship  is  initially  at  a 
speed  of  about  30  knots  with  full  pitch  ratio  (1.6).  The 

*Coastdown  as  used  hereafter  implies  setting  fuel  rate  to  idle. 
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6.  RESULTS  (Cont'd) 

6*6     Crashback  from  30  knots  (LH2500-A)  (Confd) 

first  simulation  that  was  tried  used  the  following  sequence 
of  events: 

a.  reduce  fuel  to  idle 

b.  decrease  pitch  to  -1«2 

c.  increase  fuel  to  9900  Ib/hr 

In  performing  this  simulation,  it  was  observed  that  the  pitch 
changing  operation  took  a  rather  lengthy  time.  The  end  result 
would  have  been  a  time  to  stop  in  excess  of  60  seconds.  The 
simulation  was  stopped  and  a  new  sequence  of  events  was 
attempted, 

a.  simultaneously  decrease  fuel  to  idle  and  decrease 
pitch 

b.  when  pitch  equals  -.8,  increase  fuel  to  9900 
Using  this  sequence  of  events, 'a  stopping  time  of  ^4.3 
seconds  resulted.  Results  of  this  simulation  are  plotted 
in  figures  (6~2&N  -  (6=2*0.  A  copy  of  the  computer  program 
used  in  sequence  (a)  is  included  in  Appendix  G. 

6*7     Crashback  from  9  knots  (KV  Major  12) 

The  same  sequence  of  events  described  in  section  6Ji 
were  used  here,  viz., 

1.  coastdown  -  opening  clutch 

2.  engaging  reverse  clutch 

3.  stopping  the  ship 

Forward  clutch  disengagement  occurred  at  11  seconds,  and 
6o2  seconds  later,  the  reverse  clutch  engaged,  The  ship 


51 
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6.7  Crashback  from  9  knots  (KV  Major  12)  (Cont'd) 

Was  stopped  in  95  »2  seconds  with  a  head  reach  of  1208.7  feet. 
Results  of  this  computer  simulation  are  plotted  in  figures 
(6-£.5)   -  (6-17) »     Copies  of  the  programs  used  in  all  three 
sequences  are  included  in  Appendix  H« 

6.8  Crashback  from  12  knots  (B  A  W  7K98FF) 

This  simulation  successfully  illustrated  the  increase 
in  prime  mover  speed  as  the  CRPP  went  through  a  zero  torque 
coefficient  region.  The  sequence  of  events  employed  here  were: 

1.  simultaneously  decrease  fuel  to  idle  and  pitch  ratio 

2e  when  pitch  ratio  equaled  -.2,  increase  fuel  to  stop 
the  ship 

The  initial  simulation  started  from  a  pitch  ratio  of  1.0  and 
decreased  pitch  to  -1.2.  At  this  pitch  ratio,  the  propeller 
torques  were  so  high  that  the  equivalent  of  diesel  stalling 
occurred.  The  change  in  pitch  ratio  was  modified  to  go  from 
+1.  to  -1.  This  proved  to  be  satisfactory.  The  time  of  pitch 
change  was  3^»5  seconds.  Ship  stopping  time  was  90*5  seconds 
with  a  head  reach  629*8  ft.  The  results  of  this  simulation 
are  plotted  in  figures  (6-28)  -  (6-31).  A  copy  of  the  pro- 
gram used  in  sequence  (2)  is  included  in  Appendix  I. 
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7.    CONCLUSIONS 

Simple  models  were  found  to  exist  that  enable  one  to  model 
ship  propulsion  drive  trains.  Some  of  the  models  had  to  be  refined 
by  the  addition  of  time  lags,  delays,  etc.  The  combination  of 
drive  train  components  must  be  compatible  with  the  characteristics 
■  of  each  drive  train  member.  The  simulations  illustrated  the 
importance  of  proper  sequencing  of  propulsion  plant  control  para- 
meters. Conditions  of  overspeeding,  excess  torques,  and  high 
energy  dissipation  rates  in  clutches  can  be  easily  simulated. 
These  simulations  can  prove  to  b3  immensely  useful  to  the  control 
engineer  in  helping  him  synthesize  a  controller  that  will  keep  a 
system  within  its  region  of  safe  operation. 
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8«    RECOMMENDATIONS 

1.  A  study  should  be  taken  to  determine  what  is  the  optimal 
fuel  flow  (setting)  -  pitch  ratio  combination  for  a  particular 
prime  mover  and  the  CRPP  presented  here. 

2.  A  model  of  reversing  the  B  &  W  7K93FF  should  be  developed 
and  a  crashback  simulation  of  the  7K98FF  with  a  fixed  pitch 
propeller  should  be  performed. 

3*  A  model  controller  to  keep  shaft  rpm  nearly  constant  while 

changing  pitch  ratio  on  the  CRPP  could  be  synthesized  for 

each  one  of  the  prime  movers. 
^.  A  torque  converter  model  should  be  developed  to  be  used  with 

the  prime  movers  developed  here,  and  maneuvers  should  be 

simulated. 
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APPENDIX  A 
Gas  Turbine  Models 

1.         Lurye  Model 

The  following  model  of  a  marine  gas  turbine  was  proposed  by 
Lurye  in  Reference  (10)  in  the  form  of  linearized  differential 
equations: 

Tarbocomprs  s s or 

Free  Turbine 

(T2  .  +  1)^2  =  KzlA  GT  £*+   K22A"i  -  K,  A/> 


T,  =         v>T  ....  tar  Decompressor 


±!k.  +    *n« 


time  constant 


wll 


ets/t.  .  #6.  amplification  coefficient 

of/^,*"""   o*/^         of  turbocctnpi^essor  by  fual 


K-_  =      c/^g   ••••  load  coefficient  by  the  rate 


12 


of  rotation  of  free  turbine 


T   "        ^ ••••  time  constant  of  free 

2    ""  ejM.  ,jm^  turbine 


. . . e  amplification  coefficient 

21      ,  n  \m  of  free  turbine  by  fuel 
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APPENDIX  A  (Cont'd) 


1.   (Cont'd)        ^ 


Kv  =     ° r ....  load  coefficient  by  propeller 

o>AV     dM^  pitch 

-  —3  /-  ! 

c>MT 

'X. 

K22  =      JoTp^       «...  amplification  coefficient  by 

tJfl  J/V        rate   °^  comPressor  rotation 

—     Ti 


CaJi     ~   turbocompressor  rate  of  rotation 
^2.  =  free  turbine  rate  of  rotation 
G_  =  fuel  consumption 

ft    -  propeller  pitch 

Z    =  time  lag  in  combustion  chamber 

*i-    ^1  ("Ifirx)   torque  moment  of  turbocompressor 

^2.        Tz   V^*;  T)  torque  moment  of  free  turbine 

J.  =  turbocompressor  rotational  inertia 

J  -  free  turbine  rotational  inertia 

M      =  f     {°J^Jft  )  resistance  moment  of  propeller 

m      =  f     (w/  cuj)  resistance  moment  of  compressor 
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APPENDIX  A  (cont'd) 

2.    An  Extension  of  the  Saravanamuttoo  Model 

Using  the  technique  proposed  by  Saravanamuttoo  in  References 
(7)  &  (&)>   the  model  can  be  extended  to  the  case  where  2-spool  gas 
generator- free  turbine  combination: 
HP  Rotor  equation  y  fJjSS. 

**,  K,r  ^  (-T+-  zj  -    sjr   ^  <£w  (r3~Z)  *  ^   "^oF^ 

LP  Rotor  equation  /^*    • 


/ 


Free  Turbina  equation  ^^ 


?7«-WFT  C/7  FT  (T6  "  V     " 


FT 


ctX 


IP-HP  Intercoropressor  Volume 

W        -  W        =    RT2      d '  % 
WHPC       WIPC        ~      -7-7 

Combustor[>  volume 


WHPT  *  WHPC  -  W# 


Interturbine  Volume  (IP-HP) 
Interturbine  Volume  (HP-FT) 


WHPT    FT 


k       </*" 


APPENDIX  A  (Cont'd) 


2.    (Cont'd) 
Combustor 


T^  =  usn  Nl/^C^  (V„c  5" 


C?c   (*&   flSA,r) 


Nomenclature 


W  =  mass  flow  rate 
T  =  temperature 

C  =  specific  heat  at  constant  pressure 
?Lr-   turbine  over-all  efficiency 
"%c  =  compressor  over-all  efficiency 
K  =  gas  constant 
HV=  fuel  heating  value 
V  =  volume 
N  =  rotational  speed 
P  =  pressure 
t  =  time 

I  =  polar  moment  of  inertia 
Station  Numbering 

1  =  LP  Compressor  inlet 

2  =  HP  Compressor  inlet 

3  =  corabustor  inlet 

**  =  HP  turbine  inlet  (combustor  cutlet) 

5  =  LP  turbine  inlet 

6  =  FT  inlet 

7  =  FT  outlet 
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APPENDIX  A  (Cont'd) 

2.    (Cont'd) 

Subscripts 

LPC  =  lew  pressure  compressor 
HPC  -  high  pressure  compressor 
FB  =  Corcbustor 
KPT  =  high  pressure  turbine 
LPT  =  low  pressure  turbine 
FT  =  free  turbine 
HP  =  HPC-HPT  rotor 
LP  =  LPC-LPT  rotor 


9'+ 
APPENDIX  A  (Cont'd) 

FT^A-2  &  LM2500-A  Marine  Gas  Turbines 

Figure  (A-l)  is  the  power-speed  curve  for  constant  SFC  values. 
The  torque-speed  characteristics  shown  in  figure  (A-2)  can  be 
derived  from  this  curve.  Figure  (A-3)  and  (A-*0  are  the  corres- 
ponding curves  for  the  LM2500-A.  Figures  (A»2)  and  (A-*0  show 
examples  of  the  linear  approximations  for  constant  fuel  rate 
values  used  to  approximate  the  gas  turbines*  characteristics.  For 
intermediate  values  of  fuel  rate,  linear  interpolation  was  used* 
The  subroutines  used  for  the  FT^A-2  and  the  LM2500-A  are  included 
in  the  rear  of  this  appendix, 

A  sketch  of  the  FT^A-2  with  a  JFC  25-28  fuel  controller  is 
shown  in  Figure  (A-5)<=  A  controller  similar  to  the  JFC  25-28  is 
the  JFC  25-7  shown  in  figure  (A~6).  This  controller  is  typical 
of  the  non-linear  systems  engineers  are  often  asked  to  model. 
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SUBR 

THIS 

RPM 

MINI 

MAXI 

E  N  T  E 

AND 

FT-L 

SEAR 

I  F  (  W 

I  F  (  w 

IF(W 

I  F  (  W 

I  F  (  W 

I  F  (  W 

I F  ( W 

I  F  (  W 

I F  ( W 

T0R0 

RETU 

CALL 

CALL 

TORQ 

RETU 

CALL 

CALL 

TORQ 

RETU 

CALL 

CALL 

TORQ 

RETU 

CALL 

CALL 

TORQ 

RETU 

CALL 

CALL 

TORQ 

RETU 

CALL 

CALL 

TORQ 

RETU 

CALL 

CALL 

TORQ 

RETU 

CALL 

CALL 

TORQ 

RETU 

END 

SUBR 

A00  = 

RETU 


OUTINE 

SUBROU 
AND  FUE 

run  fue 

MUM  FUE 
RING  AR 
POWER  T 
BS 


LM2  50  (RPM,WF,T0RQ,A08,A16,A2,A3,A4,A5,A6, A7,A99) 

TINE  CALCULATES  FREE  TURBINE  TORQUE  AS  A  FUNCTION  OF 

L  FLOW 

L  <ATE  IS  BOO  LP/HR 

L  RATE  ALLOWABLE  IS  9900  LB/HR 

GUMENTS  FOR  THIS  SUBROUTINE  ARE  FUEL  RATE  (WF) 

URBINE  RPM  (RPM)   OUTPUT  IS  POWER  TURBINE  TORQUE  IN 


CH 
F 

F 
F  , 

F  , 
F 

F  , 
F  , 
F  , 
F  , 


TO  B 
FQ. 
LE. 
LE. 
LF. 
LE. 
LE. 
LE.6 
LE. 
LE. 
=-1.965 
RN 
LWF08 
LWF16 
=A08+( W 
RN 
LWF2  ( 
L  W  F  1  6 
=  A  1  6  +  (  W 
RN 
LWF3  ( 
LWF2  ( 
=A2+( WF 
RN 
LWF4  ( 
LWF3  ( 
=A3+(WF 
RN 
LWF5  ( 
LWF4  ( 
=A4+(WF 
RN 
LWF6  ( 
LWF5  ( 
=A5+(WF 
RN 
LWF7  ( 
LWF6  ( 
=A6+(WF 
RN 
L  W  F  9  9 
LWF7  ( 
=A7+(WF 
RN 

OUTINE 
-1.965* 

RN 


RACKE 
30  0.  ) 
1600. 
2  0  00. 
3000. 
4000. 
5000. 
000.  ) 
70  00. 
9900. 
*RPM  + 


T  WF 
GOT 
)  GO 
)  GO 
)GOT 

)Gcr 

)  GO 
GOTO 
)GOT 

)GOT 
5500 


AND 
015 
T016 
T017 
018 
010 
TO  1 
12 
013 
014 


THEN  PROCEED  TO  FUEL  SUBROUTINES 


] 


(RP!',A08) 
(RPM, A  16) 
F-8  0  0.  )*.00125*( A16-A08) 

RPM,A2) 
(RPM,A16) 

F-1600. )*.002  5*( A2-A16) 

RPM, A3) 
RPM,A2) 

-2000. )« .001*( A3-A2) 

RPtA,  A  4  ) 
R  P  M  ,  A  3  ) 
-3000.  )*. 001» { A4-A3) 

RPM,A5) 
R  P  M  ,  A  4  ) 
-4000. )».001*( A5-A4) 

RPM,A6) 
R  P  M  ,  A  5  ) 
-5000. )«.001*( A6-A5) 

RPM,A7) 
RPM, A6) 
-6000. )*.001*( A7-A6) 

(RPM, A99) 

RPM,A7) 

-7  000. ) *. 000  3448* ( A99-A7 ) 


LWF08  (RPM,A0'8) 
RPM+5500. 
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END 

SUBROUTINE  LWF16  (KPM,A16) 

A16=-2.94*RPM+12500. 

RETURN 

END 

SUBROUTINE  LWF2  (RPM,A2) 

A2=-3. 48* RPM +16 150. 

RETURN 

END 

SUBROUTINE  LWF3  (RPMfA3) 

A3=-4.77*RPM+23250. 

RETURN 

END 

SUBROUTINE  LWF4  (RPM,A4) 

A4=-5.933«RPM+30800. 

RETURN 

END 

SUBROUTINE  LWF5  (RPM,A5) 

A5--7. 133*RPM+39500. 

RETURN 

END 

SUBROUTINE  LWF6  (RPM,A6) 

A6=-7.5*RPM+47500. 

RETURN 

END 

SUBROUTINE  LWF7  (RPM,A7) 

A7=-8.*RPM+55300. 

RETURN 

END 

SUBROUTINE  LWF99  (RPM,A99) 

A99=-10. L5«RPM+73600. 

RETURN 

END 
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APPENDIX  B 
Diesel  Particulars 

1.    Burmeister  and  Wain  7K98FF 

Test-bed  results  of  the  B  &  W  7K98FF  are  shown  in  Figure 
(B-l).  From  this  figure,  it  was  possible  to  derive  the  ideal 
diesel  power-speed  curves  for  constant  fuel  settings;  this  is 
shown  in  figure  (B-2).  The  torque-speed  characteristics,  figure 
(B-3),  are  in  turn  derived  from  figure  (B-2).  The  "FF"  at  the 
end  of  the  designation  of  this  engine  denotes  that  this  diesel  is 
intended  for  direct  drive,  i.e.,  no  reduction  gears.  Other 
technical  data  is  presented  below: 

continuous  service  rating  2^500  metric  BMP  @  100  rpm 
maximum  continuous  service  rating  26600  metric  BHP  ©  103  rpm 
Length  *J-5»  2" 
Breadth  l6«  2" 
Height  k>y   2" 

The  computer  program  used  to  model  the  diesel  is  included 
at  the  end  of  this  appendix. 
2.    Mirrless  KV  Major  12 

A  cutaway  drawing  and  technical  specifications  of  this  engine 
is  shown  in  figure  (B-^).  The  "12"  at  the  end  of  the  designation 
indicates  the  number  of  cylinders.  Figure  (B-5)  shows  the  fuel 
map  for  this  engine.  Using  this  fuel  map  and  the  elementary 
diesel  equation  discussed  in  section  2,  it  was  possible  to  obtain 
the  approximately  Linear  torque-fuel  setting  characteristics  shown 
in  figure  (B-6).  The  relationship  obtained  from  figure  (B-6)  is 


108 


APPENDIX  B  (Cont'd) 

2.    (Cont'd) 

T  =  550z 

z  =  fuel  sotting  in  percent 

T  =  torque  in  ft-lbs 
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APPENDIX  C 
Clutches 

Figures  (C-l)  and  (C-2),  which  were  provided  by  the  Philadelphia 

Gear  Company,  describe  the  operation  of  quill  shaft  version  of  the  firn^s 

synchroclutch.  Additional  useful  information  provided  by  Philadelphia 

Gear  Is: 

Inertia Input  -17^  lb-ft* 

Output  -110  lb-jNT 

torsional  stiffness* •• »•  Approx.  270  million  in-lb/Rad 

Figures  (C-3)  and  (C-^*)*  which  were  provided  by  the  Koppers  Co., 

describe  the  operation  of  the  BIH-Dynetic  clutches. 


*  For  information  purposes  only 


U',. .:  :;  zv.ii':^ii 


117 


H        C      E        I  J 


QUILL  SHAFT 


FRICTION  CLUTCH 
COOLINC  OIL  PORTS 


A.   FRICTION  CLUTCH  PACK 

B    ROTARY  UNION 

C    DENTAL  CLUTCH  ASSEMBLY 

D.  ENGAGING  CYLINDER  FOR  DENIAL 
CLUTCH 

E.  THRUST  BRG. 

F.  SHIFT  COLLAR 

G.  PINS 
H.  SLEEVE 

I.  INTERMEDIATE  GEAR 

J.  DRIVE  SPLINE 

K.  OUTPUT  FLANGE 

I  SELF-LOCK  PISTON  ASSY. 

M.  DISENGAGE  CYLINDER-DENTAL 
CLUTCH 

N.   SPACE  FOR  OPTION.  QUILL  SHAH 
BRAKE 

0.   FRICTION  APPLY  CYLINDER 

P.   HYDRAULIC  BALANCE  &  FRICTION 
RETURN  PISTON 

Q.   JOURNAL  S  THRUST  BEARING 
ASSEMBLY 


FIG.  2    CROSS  SECTION  OF  TYPICAL  SYHCHROCLUTCN  (QUILL  SHAFT  VERSION) 


[-\:r,r:,.vr\:i\:,i.]  ivcm  rv!^,;/-,;j 


■  The  new  Philadelphia  SYNCHROCLUTCH  consists  of  a 
friction  clutch  in  parallel  with  a  positive  drive  dental  clutch. 
The  friction  clutch  synchronizes  the  driving  and  driven 
shafts,  after  which  the  dental  clutch  is  engaged  to  provide 
a  positive  drive.  The  dental  clutch  features  fully  crowned, 
nitrided  gears  for  maximum  capacity  for  a  given  diameter. 
The  friction  clutch  is  of  the  multi-disc  type  with  positive 
oil  circulation  for  long  life  and  heat  removal.  The  friction 
clutch  has  ample  capacity  to  absorb  thermal  loads  imposed 
by  the  synchronizing  cycle — even  very  high  loads  encoun- 
tered when  reversing  large  marine  gears  in  the  40,000 
plus  HP  range.  In  addition,  dead  shaft  pickup  is  possible 
in  large  marine  propulsion  systems.  This  permits  the  prime 
mover  to  be  run  at  speed  in  port  for  auxiliaries;  such  as 
generators,  while  retaining  the  ability  to  get  underway 
instantly  by  clutching  in  the  propeller. 

The  dental  clutch  features  a  unique  helical  gear  element 
which  automatically  aligns  the  teeth  of  the  dental  clutch 
prior  to  engagement.  This  feature  is  illustrated  in  Fig.  3. 
The  clutch  may  bo  actuated  hydraulically  or  pneumatically, 
using  low  pressure  air  (100  psi)  or  oil  (150  psj).  The 
synchronizing  action  is  controlled  by  a  transistorized  dif- 
ferential speed  control  using  non-contacting  magnotic 
pick-ups  for  speed  indication  and  position  readout.  Once 
engaged,  the  dental  clutch  does  not  require  power  to 
remain  in  mesh,  thus  eliminating  dependence  on  electrical 
or   hydraulic   power  for  continued  torque   transmission. 


External  jacking  is  provided  for  emergency  manual  engage- 
ment and  disengagement  of  the  dental  clutch. 

Referring  to  Fig.  2,  the  friction  clutch  (a)  is  first  engaged 
to  force  the  driving  and  driven  shafts  into  synchronism. 
The  clutch  is  clomped  by  admitting  oil  (or  air)  to  the  piston 
through  the  rotary  union  (b).  After  the  driving  and  driven 
shafts  have  reached  a  speed  differential  of  approximately 
]/2%.  the  transistorized  speed  control  signals  for  dental 
clutch  (c),  engagement  by  admitting  oil  (or  air)  pressure 
to  the  engaging  side  of  the  dental  clutch  piston  (d).  The 
piston  transmits  the  thrust  through  a  heavy  duty  thrust 
bearing  (e),  to  the  rotating  shift  collar  (f),  which  is  con- 
nected through  large  pins  (g),  to  the  sleeve  (h),  of  the 
dentai  clutch.  The  sleeve  carries  with  it  the  intermediate 
gear  (i),  and  meshes  with  the  drive  spline  ()),  on  the  output 
flange  (k),  in  the  manner  illustrated  in  Fig.  2.  A  self-locking 
piston  (I)  on  the  dental  clutch  prevents  piston  back-off 
under  torque  reversal.  This  action  is  shown  in  Fig.  3.  To 
disengage,  pressure  is  admitted  to  the  reverse  side  (m)  of 
the  dental  clutch  piston,  which  pulls  the  dental  elements 
out  of  engagement  as  shown  in  Fig.  3.  The  friction  clutch 
may  then  be  de-energized  to  permit  the  driving  and  driven 
shafts  to  rotate  relative  to  each  other.  Where  desired,  a 
quill  shaft  brake  (n)  may  be  incorporated  into  the  clutch 
to  stop  the  quill  shaft  from  rotating,  due  to  residual  drag 
torque  within  the  friction  clutch.  This  feature  requires  no 
additional  length. 


Fig*   (C-l) 


?■  8 


■ 


1.   DENTAL  CLUTCH  DISENGAGED 

INTERMEDIATE  GEAR  /DRIVE   SPLINE 


HELICAI    GEAR 


2.   INTERMEDIATE  GEAR  &  DRIVE  SPLINE  GEAR 
BUTTING 


Mm 


3.  SLIDING  SLEEVE  HELICAL  GEAR  ROTATES 
INTERMEDIATE  GEAR  INTO  ALIGNMENT  WITH 
DRIVE  SPLINE  TEETH 


1 


4.  SLIDING  SLEEVE  PUSHES  INTERMEDIATE 
GEAR  &  DRIVE  SPLINE  TEETH  INTO  FULL 
MESH. 


wrCw I  ■" I  wn  i  I  U  111  w 


JO. 

11. 


12. 


SPEED  DIFFERENTIAL  (INPUT/OUTPUT). 
ACTUATION  FiuiD. 


UP  TO  100%  OF  RATED  SPEED. 

OIL  @  150  PSI.  (SPECIAL-AIR  @  100  PSI.) 


FRICTION  MATERIAL. 


PGC.  FS001  GRADE  RC1. 


STEEL  DISC  MATERIAL. 


HARDENED  STEEL.  HIGH  CAREON.  GROUND. 


ROTATING  PARTS. 


FORGED  STEEL,  HEAT  TREATED  AND 
DYNAMICALLY  BALANCED. 


DENTAL  CLUTCH  MEMBERS. 


NITRIDED  STEEL. 


DENTAL  CLUTCH  TEETH. 


CROWNED.  CHAMFERED  ON  ENTRY  SIDE. 


ROTARY  UNION. 


EXTERNAL— UNIT  CONSTRUCTION. 


ACCESSIBILITY. 


J3JRCCTION  OF  ROTATION. 

PERMISSIBLE  AXIAL  MOTION  OF  OUTPUT 
FLANGE  CONNECTION. 


ALL  CLUTCH  MEMBERS  REMOVABLE  WITH- 
OUT  REMOVING  HUB  OR  QUILL  SHAFT. 


CW  OR  CCW. 


±1/8"  STANDARD. 


PERMISSIBLE  RADIAL  RUNOUT  AT  OUTPUT 
FLANGE. 


.015  TIR  ON  500  SIZE  AND  LARGER. 
.010  TIR  ON  120  SIZE. 


ELECTRICAL  POWER  REQUIREMENTS  FOR 
CONTROLS. 


115V.  AC-50/60  HZ,  50  WATTS. 


14. 
15. 


CLUTCH  ENGAGEABI.E  AT  STANDSTILL. 


YES. 


EMERGENCY  MANUAL  CLUTCH  OPERATION. 


YES. 


OPERATION   AND 


UNIFICATION   OF  SYNCHROCLUTCH 
Fig.    (C-2) 
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DISENGAGED  POSITION 


ENGAGED  POSITION 


FUNDAMENTAL  DESIGN 

All  BLH-Dynetic  Clutches  are  of  the  same  basic 
design,  modified  as  necessary  to  meet  a  wide  vari- 
ety of  application  conditions. 

The  forces  developed  during  engagement  are 
completely  self-contained,  and  produce  no  thrust 
loads  on  the  shafts  or  clutch  housing.  After  engage- 
ment, power  is  transmitted  through  hardened  and 
nitrided  coupling  teeth  in  positive  mesh,  with  no 
loss  due  to  slippage.  Hydraulic  pressure  is  required 
only  during  the  actual  engaging  or  disengaging 
cycle.  The  parts  are  held  in  their  final  positions  by 
spring  loaded  detents,  or  mechanical  locks,  as 
necessary. 


MAJOR   COMPONENTS 

]  Operating  Cylinder— to  provide  the  axial  force 
necessary  to  engage  or  disengage  the  clutch.  Actu- 
ation is  usually  hydraulic  using  lubricating  oil  as  the 
fluid,  but  pneumatic  operation  is  also  available. 
[..^  Friction  Discs— to  develop  friction  torque  neces- 
sary to  synchronize  the  shaft  speeds.  These  discs 
are  cooled  by  a  continuous  flow  of  lubricating  oil. 
[2j  Transition  Torque  Control— spring  loaded  pins 
key  with  lugs  in  the  load  end  hub  to  transmit  torque 
during  the  transition  from  friction  to  gear  tooth 
drive,  and  to  align  gear  teeth  prior  to  engagement. 

Q  Gear  Coupling— to  transmit  power  after  engage- 
ment is  complete  and  accommodate  shaft  misalign- 
ment, if  present. 


DESCRIPTION  OF  BLII  DYNETIC  CL'JTC:* 


Fig.  (C-3) 
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Figure  A 


10.10 


• 


SP  V 


Figure  B 


(   :H 


SP 


Figure  C 


SP 


n 

Figure  D 


The  heart  of  the  synchronizing  mechanism  is  the 
cluster  of  synchronizing  pins  (SP)  which  cooperate 
with  the  holes  in  the  power  end  hub  (PEH)  through 
which  they  project  to  control  engagement  as  indi- 
cated by  the  diagrams  at  the  left. 

As  engagement  begins  the  actuating  piston  moves 
the  power  end  hub  and  its  connected  parts  toward 
the  load  end  hub.  Through  detent  action,  an  initial 
friction  torque  is  developed  which  swings  the  small 
diameters  of  the  synchronizing  pins  against  the 
sides  of  their  respective  holes  in  the  power  end  hub, 
as  in  Figure  A. 

Continued  pressure  brings  the  cam  surfaces  of 
the  synchronizing  pins  into  contact  with  those  of 
the  power  end  hub,  thereby  transmitting  full  piston 
pressure  to  the  friction  disc  pack  as  in  Figure  B. 
The  friction  torque  so  developed  acts  to  synchronize 
shaft  speeds.  The  geometry  of  the  assembly  is  such 
that  the  parts  remain  in  the  position  shown  as  long 
as  slipping  occurs  in  the  friction  discs. 

Slipping  stops  at  the  friction  discs  when  the  shaft 
speeds  synchronize.  This  unbalances  the  forces  and 
permits  the  power  end  hub  to  advance  over  the  syn- 
chronizing pins  as  in  Figure  C.  At  the  same  time, 
spring  loaded  pins  are  moved  into  spaces  between 
lugs  in  the  end  of  the  load  end  hub,  acting  as  loose 
keys  to  control  torque  until  the  coupling  teeth 
engage.    ■ 

As  the  synchronizing  pins  become  aligned  with 
their  mating  holes  as  in  Figure  D,  the  cam  surfaces 
separate,  removing  the  axial  force  on  the  friction 
disc  pack.  This  permits  the  spring  loaded  pins  to 
contact  the  sides  of  their  lugs,  maintaining  syn- 
chronism of  shaft  speeds  and  aligning  the  coupling 
teeth  for  engagement  as  the  piston  completes  its 
travel. 

With  completion  of  engagement  the  hydraulic 
pressure  is  cut  off,  and  the  assembly  transmits 
power  as  a  simple  gear-type  coupling. 

During  disengagement  the  piston  moves  the  parts 
away  from  the  load  end  hub,  disengaging  the  spring- 
loaded  pins  from  their  lugs  before  the  coupling 
teeth  disengage,  and  securing  the  assembly  in  the 
final  position  to  prevent  dynamic  unbalance. 


OPERATION    OF    13LH   DYNETIC   CLUTCH 


Fig.    (C-*0 
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APPENDIX  D 
Propeller  Subroutines 

This  appendix  contains  the  computer  programs  used  to  model  the 
B*4--?0  (p/d  =  1.)  propeller  and  Read's  CRPP,  both  of  which  are  discussed 
in  Section  ^.  Included  in  this  appendix  is  the  program  described  in 
Section  5»^  that  computes  shaft  friction* 
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SUBROUTINE  CRPP  (  J  ,  C0»  CT  ,  AP  1  ,  AP8  ,  BP  1 ,  BP8  ,  AMO  ,  AMI ,  BMfi  ,  RM1 ,  AM4  ,  BM'i ,  A 
2M2  t  BM2  i  AMO  ,  BMO  ,  AP't ,  B P4  ,  PD  ) 
CRPP  COMPUTES  MODIFIED  THRUST  AND  TORQUE  COEFFICIENTS  FOR  A.  GIVEN 
VALUE  OF  MODIFIED  ADVANCE  COEFFICIENT  AND  PITCH  DIAMETER  RATIO. 

PROPELLOR  DATA  IS  TAKEN  FROM   DYNAMIC  ANALYSIS  AND  SIMULATION 
OF  SHIP  AND  PROPULSION  PLANT  MANUVERING  PERFORMANCE   BY  WENDEL 
AND  DUNNE  APPEARING  IN  PROCEEDINGS  OF  THE  SECOND  SHIP  CONTROL 
SYSTEMS  SYMPOSIUM. 


MODIFIED  ADVANCE  COEFF J=VA/SQRT( VA»*2+ ( N«D ) **? ) 

MODIFIED  THRUST  COEFF. . .CT=T/(RHO«(D*»2)*(VA**2+(N*D) *»2) ) 
MODIFIED  TORQUE  COEFF. ..CQ=Q/(RH0*(D**3)*(VA*»2+(N*D)**2)) 


RHO  =  DENSITY  OF  WA'l  ER 
D=PROPELLOR  DIAMETER 
VA=PROPELLOR  SPEED  OF  ADVANCE 
N=PROPELLOR  ROTATION  RATE 
T=THRUST 
0=TOROUE 


PITCH  DIAMTER  RATIOS  ARE  LIMITED  BETWEEN  +1.6  TO  -1.2 
ENTERING  ARGUMENTS  ARE  MODIFIED  ADVANCE  COEFFICIENT   (J) 
PITCH  DIAMTEER  RATIO  (PD) 


AND 
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OUT 
REA 
IF( 
IF( 
IF( 
IF( 
IF( 
IF( 
CAL 
CAL 
F=( 
C0  = 
CT  = 
RET 
CAL 
CAL 
F=( 
CQ  = 
CT  = 
RET 
CAL 
CAL 
F=( 
C0= 
CT  = 
RET 
CAL 
CAL 


PUTS 
L  J 


ARE  CO  AND  CT 


LE 
LE 
LE 
LE 
LE 
LE 
^L 
PL 


J-.8 

AP8  + 
BP8  + 
URN 
L  MI 
L  MI 
J  +  l. 
AM1  + 
BM1  + 
URN 
L  Ml 
L  MI 
J+.8 
A  M  8  + 
8  M  8  + 
URN 
L  MI 
L  MI 


.-.8)GOT010 
.-.'*)  GOTO  11 
)G0T012 

G0T013 
G0T014 
JG0T015 
(PD,APltBPl) 

(  P  D  ,  A  P  8  ,  B  P  8  ) 

P1-AP8) 
P1-HP8  ) 

( PDtAMl,BMl) 

PD,AM8,BM8) 


.-.2 
.0.  ) 
..4) 
.  .8 
US1 
US8 
)*5, 
F»(  A 
F  *  (  B 

ONE 
N8  ( 
)*5. 
F*(  A 
E*(B 

N4  ( 
N8  ( 
)*5. 
F*(A 
F  *  (  B 


N2 

NA 


M8-AM1 ) 

M8-BM1 ) 

PD, AM4,BM4) 
PD,AM8,BM8) 

M4-AMR ) 
MA-BM8 ) 

PD,  AM2,BM2) 
PD,AM4,BMA) 
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15 


10 
11 


10 
11 


12 


13 


10 

11 
12 


F=(J+.4)« 

CQ=AM4+F* 

CT=BM4+F* 

RETURN 

CALL  ZERO 

CALL  MIN2 

F=( J+.2) * 

CQ=AM2+F« 

CT=BM2+F« 

RETURN 

CALL  PLUS 

CALL  ZERO 

F=J*2.5 

CO=AMO+F* 

CT=BMO+F* 

RETURN 

CALL  PLUS 

CALL  PLUS 

F=( J-.4)» 

CQ=AP4+F* 

CT=3P4+F* 

RETURN 

END 

SUBROUTIN 

I  F  (  P  0  .  G  E  . 

APl=-.477 

GOTOU 

API =.00 51 

BPl=-.4 

RETURN 

END 

SUBROUTIN 

I  F  (  P  D .  G  E . 

AP8=-.06* 

G0T011 

AP8=.026* 

IF(PD.LT. 

IF(PD.GT. 

BP8=.0R5« 

RETURN 

3PR  =  .  165* 

RETURN 

BP8=.243* 

RETURN 

END 

SUBROUTIN 

IF(PD.GT. 

IFCPD.LT. 

AP4=.015 

GOTO  12 

AP4=.0843 

GOTO  12 

AP4=-.123 

IF(PO.LT. 

IF(PD  ,GT 

F'P4  =  .24  5» 


( AM2-AM4) 
(BM2-BM4) 

( PDtAM0,BM0) 
{  PDf  AM2f.BM2) 

5 . 

(AM0-AM2 ) 

(BM0-BM2 ) 

4  (PD,AP4,BP4) 
(PDf AMO.BMO) 

(AP4-AM0) 
(BP4-8K0) 

3  (PD,AP8,BP8) 

4  ( PD, AP4,BP4) 
2.5 

(AP8-AP4) 
(BP8-BP4) 


E  PLUS1  (PD,AP1,BP1) 

.63JG0T010 
*PD-.015 

6*PD-.049 


E  PLUS8  (PD,AP8tBP8) 
. 35)GOT010 

PD 

PD-.03 

-.A) GOTO  12 

.4)  G0T013 
PD-.275 

P.D-.243 

PD-.3382 


E  PLUS4  (PD,AP4,RP4) 
.4) GOTO  10 
-.4) GOTO  11 


*PD-.019 

*PD-. 032 5 

-.4)G0T013 
.  .4)  G0T014 

PD-.9 
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11 
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RETURN 
BP4=.37 
RETURN 
BP4=.41 
RETURN 
END 

SUBROUT 
IF( PD.L 
I F  (  P  D  .  G 
AM0=.01 
G0T012 
AM0=-.  1 
GOTO 12 
AM0=. 10 
I  F  (  P  D  •  L 
1F(PD.G 


BMO  = 


1 


RETURN 

3M0=.53 

R  E  T  U  R  N 

BM0=.53 

RETURN 

END 

SUBROUT 

IF(PD  . 

I  F  (  P  D  . 

AM2=.01 

G0T012 

AM2=.  12 

G0T012 

AM2=-.l 

G0T012 

IF(PD  . 

IF(PD  . 

BM2=.31 

RETURN 

BM2=.53 

RETURN 

BM2=.53 

RETURN 

END 

SUBROUT 

IF(PD.L 

IF(PD  . 

AM4=.00 

G0T012 

AM4=-.0 

G0T012 

AM4=.07 

I  F  (  P  0  . 

I  F  (  P  D  . 

BM4=.  18 

RETURN 

BM4=.5« 

RETURN 

BM4=.4* 


*PD-.04 
5*PD-. 150 


INE  ZERO  (PD,AMO,BMO) 

T.-.4JGOT010 

T.  .4) GOTO  11 


25*PD-.0375 

92*PD-.031  1 
T.-.625)C0T013 
T.  .425)G0T014 

0*PD+.020 

*P0+. 1575 

* PD-. 07 32 

INE  MIN2  (PD,AM2,BM2) 

GT.  .045) GOTO  10 

LT.  -.0475)GOT011 

25 

3*PD-.043 

32*PD-.05 

LT.  -.625)G0T013 
GT.  .425JG0T014 
0*PD+.02 

*PD+.1575 

•PD-.0732 


INE  MIN4  (  PD,  AM4tBy,4  ) 

T.-.425)GDT010 

GT.  .275)GOT011 

6 

925*PD~.0335 

25*PD-.014 
GT.  .88)  G0T013 
LT.  -.36J0OTO14 
4*PD+.076 

PD-.202 

PD  +  .  145 
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RETURN 

END 

SUBROUTINE  MIN8  (PD,AM8,BMR) 

IF(PD  .LT.  -.55) GOTO  10 

AM8=.07*PD 

GOTO  13 

AM8=-.028*PD-.054 

IFtPD  .LT.  -.45)GOT011 

42)  G0T01? 

2760 


I  F  (  P  D  .  G  T  . 

BM8=.09*PD+ 

RETURN 

11  BM8=.21*PD+.33 
RETURN 

12  BM8=.283»PD+. 1652 
RETURN 

END 

SUBROUTINE  MI  ONE  (PD,AM1,BM1) 

AM1=.064*PD-.018 

IF(PD  ,GT.  .6)G0T011 

IF(PD  .LE.  .400) GOTO  10 

BM1=-.88*P0+. 743 

RETURN 

10  BM1=.391 
RETURN 

11  BMl=.215 
RETURN 
END 


123 


PAGE   1 


05/03/ 


SUBROUTINE  SHAFT  (SN,QF) 

THIS  SUBROUTINE  COMPUTES  SHAFT  FRICTION 

ENTERING  ARGUMENT  IS  SHAFT  ROTATION  RATE  (SN)  IN  REV  PER  SECOND 
OUTPUT  IS  FRICTION  TORQUE  (QF)  IN  FT-LBS 

SN=ABS( SN) 

IF(SN.LE. .417JGOT010 
0F=SN*6000. 
RETURN 
10  QF=2500. 
RETURN 
END 
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APPENDIX  E 
Ship  Dynamics 

1,  Ship  Propulsion  Equations 

The  ship  propulsion  equations  with  no  external  forces  or 
.  rudder  movement  are  the  two  first-order,  non-linear,  differential 

equations  given  below: 

civ 
m  -rr  =  T  -  R  (thrust  equation) 

dt  =  QD  "  QF  "  %       (tor(ine   equation) 

v  =  ship  velocity  (ft/ sec) 

T  ~  thrust  (lbs) 

R  =  resistance  (lbs) 

ra  =  mass  of  ship  plus  10$  (1D^C  ) 

I  -  rotational  inertia  of  drive  shaft  (Ib-ft-sec  ) 

Q  =  prime  mover  drive  torque  on  propeller  (Ib-ft) 

Q  =  shaft  friction  torque  (Ib-ft) 

o  =  propeller  torque  (lb«ft) 

n  =  propeller  angular  speed  (rev/ sec) 
The  two  equations  are  coupled  through  the  propeller  terms  T  and 
Q  •  The  equations  were  solved  simultaneously  on  an  IBM  1130 
computer  using  a  fourth  order  Runge-Kutta  integration  technique, 

2.  Sign  Convention 

The  sign  convention  used  here  is  that  propeller  angular  speed, 
n,  and  velocity,  v,  are  defined  as  positive  for  forward  ship  motion. 
Torques  are  considered  positive  when  acting  in  the  direction  of 


130 
APPENDIX  E  (Cont'd) 

2.    (Cont'd) 

positive  n;  thrusts  are  considered  positive  in  the  direction  of 

positive  v. 
3»  •  Ship  Displacement 

For  all  simulations,  a  ship  with  a  displacement  of  ^000  tons 

was  used.  In  the  equations  of  motion,  10$  of  this  figure  was 

included  to  account  for  entrained  water. 
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APPENDIX  F 
FT^A-2  Drive  Train 

1.    Digital  Simulation  of  an  FT^A°2  Forward  Drive  Train 

Figure  (F-l)  is  a  sketch  of  a  drive  train  powered  by  two 
FT4A-2  marine  gas  turbines  that  was  used  in  the  simulations. 
Except  for  the  propeller  the  drive  train  is  basically  the  same 
used  by  Rubis  in  Reference  (11).  The  reduction  gear  ratios  relate 
propeller  speed  and  torque  to  power  turbine  speed  and  torque. 

The  drive  train  moment  of  inertia  was  assuned  to  be  the  same 

<;        2 
as  Rubis S  viz.,  2.88  x  l(r  lb»ft~sec  (referred  to  propeller 

shaft). 

It  should  be  noted  that  at  high  speeds  and  low  fuel  rates,  the 
FTAA-2  gas  turbine  subroutine  will  yield  negative  torques.  A 
statement  was  made  in  the  main  program  to  ignore  this  and  sub- 
stitute instead  a  linear  windage  torque  of  the  form, 

Q  =  -2N  +  3000  (ft-lb) 

N  -  Free  Turbine  Speed  (RPM) 

For  the  crashback  simulation,  a  Falk  clutch  whose  torque 
capacity  is 

Q„  =■■  &S7.5         Pnet 
where 

*net  =  net  air  pressure  acting  on  shoes 
was  used  (cf.  section  3«2).  A  supply  air  pressure  of  200  psi  was 
assumed.  The  criteria  for  clutch  disengagement  was  the  same  used 
by  Rubis  in  Reference  (11):  whenever  the  propeller  torque  "seen" 
by  the  clutch  exceeded  clutch  torque,  clutch  disengagement  was 
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1.    (Cont'd) 

was  assumed  to  have  taken  place  instantaneously.  For  the  crash- 
back  simulation,  the  equation  for  forward  clutch  disengagement 
torque  was 

Q„  =  56S7.5   (195-30WW2)  (ft-lb) 

t  =  time  (sec) 

N  =  propeller  shaft  rotation  rate  (rev/sec) 
The  propeller  torque  opposing  the  forward  clutch  was 


^pc  "  -5— 


K3     &.3197 

K~  s  reduction  gear  ratio  (see  fig.  (?-!)) 

0   =  propeller  torque  ft-lbs 
During  reverse  clutch  engagement,  reverse  clutch  torque  is  given 
by 

QCR  =  5687.5  (5t-5~^n2)  (ft-lb) 

t    =  time  (sec) 

n    =  propeller  shaft  rotation  rate  (rev/sec) 
The  open  drive  train  is  defined  to  exist  whenever  the  quantity 
( 5t- 5-^-n  )  is  negative;  if  this  occurs,  Q„„  is  set  to  a  zero 
value  obviating  a  negative  value  for  clutch  torque.  A  slip  speed 
across  the  clutch  exists  during  reverse  clutch  engagement  and 
this  is  given  by 

Hslip  -  K  N  -  (Aj 


T*2 

N  =  propeller  rotation  rate  (rpm) 

N~  -  power  turbine  speed  (rpm) 
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APPENDIX  F  (Cont'd) 


1.    (Cont'd) 


K, ,  K2,  K~,  =  reduction  gear  ratios 


The  product  of  N  v  and  Q   is  the  power  dissipated  in  the  clutch, 

Si-lp       OK 

and  the  integral  of  power  over  time  is  a  measure  of  the  energy 
absorbed  in  the  clutch. 

Reverse  clutch  engagement  was  assumed  to  exist  the  instant 
N  ,.  equaled  zero. 


135 


APPENDIX  G 
LM2500-A  Drive  Train 

!•    Digital  Simulation  of  an  LH2500-A  Powered  Drive  Train 

Figure  (G-l)  is  a  sketch  of  the  drive  train  powered  by  two 
LM2500-A  marine  gas  turbines  with  a  controllable  reversible  pitch 
propeller  (CRPP).  A  popular  control  strategy  is  to  increase  speed 
by  first  increasing  propeller  pitch  while  keeping  shaft  rotation 
rate  constant;  once  full  pitch  is  attained,  shaft  rotation  rate 
is  increased.  (2*0  A  proportional- integral  controller  of  the  form 

_wy_  =  5  +  hi 

RPM        S 

W_  =  fuel  rate  (Ib/hr) 
v 

RPM  =  (ordered  shaft  RPM)- (Actual  Shaft  RPM)   (RPM) 

s   =  differential  operator 
was  used  in  one  simulation  to  keep  shaft  rotation  rate  constant 
while  increasing  pitch e 

As  in  the  FT^A-2  subroutine,  the  LM2500-A  subroutine  will 
also  yield  negative  torques  for  high  power  turbine  RPM  and  low 
fuel  rates.  These  torques  were  also  ignored  and  a  windage  torque 
expressed  by 

^  =  -1.965  N  +  5500  (ft-lb) 

N  =  power  turbine  speed  (RPM) 
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SUBR 
ACCE 
COMM 
VA=. 
CALL 
DENO 
COMP 
Y(5) 
COMP 
Y(6) 
PROP 
Y(6) 
ETA 
ETA  = 
Y(7) 
Y(7) 
RT  I 
RT=1 
SHAF 
SN  =  Y 
CALL 
AN3  = 
W  F  =  Y 
FT4A 
CALL 
DELA 
Q3  =  0 
Q3  =  Q 
03  =  0 
Y(9) 
Y(9) 
QD  I 
QD  =  2 
Y(  1  ) 
DY(  1 
Y(8) 
Y(8) 
Y(2) 
DY(2 
Y(3) 
DY(  3 
Y(4) 
Y(A) 
Y(  16 
Y{  16 
Y(  17 
Y(  17 
IF(T 
RETU 
END 


OUTI 
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UTE 
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SHA 
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COMM 

VA=. 

CALL 
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APPENDIX  H 
KV  Major  12  Drive  Train 

Figure  (H-l)  is  a  sketch  of  the  drive  train  used  for  simulation. 

2 
The  WR  of  the  engine-flywheel  was  provided  by  the  Hawker-Siddeley  Co. 
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SURROUTI 
COASTDOW 
COMMON  T 
Y(  1  )  IS 


Y(2) 
Y(3) 
Y(A) 

Y(5) 


IS 

IS 
IS 
IS 


VA=.96*Y 
CALL  PRP 
DENOM=VA 
COMPUTE 
Y(5)=262 
COMPUTE 
Y(4)=171 
PROPFLLO 
Y(4 )=.98 
ETA  IS  R 
ETA=.75+ 
Y  6)  IS 
Y(6)=.59 
Y(7)  IS 
Y(7)=97. 
IF(Y(7) 
TORO  COM 
FUEL  SET 
TORQ=550 
RT  IS  SH 
RT=104.8 
SHAFT  CO 
SN  =  Y( 1) 
CALL  SHA 
Y(8)  IS 
Y(8)=T0R 
OD  IS  DI 
00=3.72* 
DY ( 1 ) = ( 0 
DY(2)=( Y 
DY(3)=Y( 
OCF  IS  F 
QCF=6000 
QSLK  IS 
0SLK=Y(5 
IF( APS(0 
Y(9)=0CF 
Y(  10)=0S 
RETURN 
END 
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APPENDIX  I 
B  &  W  7K98FF  Drive  Train 

Figure  (1-1)  is  a  sketch  of  the  7K98FF  drive  train.  The  WR2 
of  the  engine  was  not  available  in  the  literature  surveyed.  An 
official  of  the  Burmeister  and  Wain  Corporation  provided  the  following 
weights : 

Piston  •••••  ^.18  tons 

Crosshead...  3*35  tons 

Crank  ^.97  tons 

From  the  "Shock  and  Vibration  Handbook",  Vol.  3,  p.  38-2,  the  following 

2 
expression  was  found  for  finding  the  WR  of  a  diesel: 


(  2   +  Wc  (1"  2"  }  ) 

W  =  weight  of  piston,  piston  pin,  and  cooling  fluid  (lb) 

W  =  weight  of  connecting  rod  (lb) 

h  =  fraction  of  rod  length  from  the  lower  center  of  gravity 

R  =  crank  radius  (in) 
From  the  above  expression,  the  polar  moment  of  inertia  was  estimated 
at  *U25  x  lCp  lb-ft-sec2. 
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